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SECTION  1 


INTRODUCTION  AND  SUMMARY 

This  study  has  sought  to  define  parametrically  the  capability  of  a 
high-flying  spectrometer  system  to  detect  and  characterize  particular 
gaseous  emissions  by  ground- level  stationary  sources.  The  systems  con¬ 
sidered  utilize  a  Fourier  Transform  Spectrometer  (FTS),  although  many  of 
the  study  results  are  applicable,  with  suitable  redefinition  of  parameters, 
to  non-multiplexing  spectrometers.  The  baseline  systems  employ  mosaic 
focal  planes  together  with  optical  designs  or  data  processing  techniques 
that  effect  "background  subtraction".  This  is  a  basic  requirement  of  the 
sensor  because  the  target  spectral  emission  or  absorption  is  but  a  small 
perturbation  in  the  normal  background  spectrum.  Detection  and  quantifica¬ 
tion  are  based  on  the  correlation  of  the  observed  target  spectral  contrast 
with  predicted  target  spectral  signatures.  Spectral  correlation  detection 
implies  a  resolution  requirement  (resolved  spectral  interval)  of  the  order 
of  0.1  cm  which  is  the  approximate  halfwidth  of  spectral  lines  in  the 
target  emissions  and  atmospheric  interferences. 

The  system  must  have  a  very  small  instantaneous  f ield-of-view  to  isolate 
the  target,  and  a  telescope  large  enough  to  collect  sufficient  power  for 
detection.  For  the  IFOV's  and  telescope  sizes  envisioned,  the  Afi  product, 
or  system  throughput,  is  small  compared  to  that  achieved  in  other,  more 
common,  FTS  applications;  i.e.,  the  system  throughput  can  be  accommodated 
by  a  Mlchelson  cube  with  dimensions  of  one  or  two  inches.  However,  the  FTS 
system  still  enjoys  a  considerable  throughput  advantage  over  a  grating 
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spectrometer  system.*  The  system  detectors  can  achieve  background- limi ted 
performance  under  typical  operating  conditions,  which  negates  the  FTS 
"multiplex  advantage". 

A  related  hardware  development  feasibility  study  performed  by  SS6,  Inc 
compared  both  FTS  and  grating  spectrometer  concepts  as  candidates  for  this 
application.  The  study  showed  that  when  the  relative  throughputs  and  effi¬ 
ciencies  and  detector  array  requirements  are  compared,  the  FTS  is  clearly 
superior,  assuming  the  detector  arrays  are  representative  of  current 
technology.  It  showed  also  that  the  FTS  system  has  sensitivity  equal  to  or 
better  than  a  grating  system  even  when  mosaic  array  technology  beyond  the 
current  state-of-the-art  is  assumed  for  the  grating  concepts. 

The  basic  radiometric  design  requirements  and  available  tradeoffs  for 

this  type  of  FTS  application  were  established  in  the  previous  "Down-Looking 

2 

Interferometer  Study".  That  study  also  evaluated  some  of  the  impacts  of 
"then-current"  technology  relative  to  detectors  and  A/D  converters,  and  dis 
cussed  the  merits  of  cooling  the  system  foreoptics.  Since  it  is  not 
generally  available,  salient  points  from  that  report  are  developed  here  as 
needed. 

The  following  were  the  specific  objectives  of  this  study: 

a.)  to  define  several  baseline  radiometric  concepts  of  a  system 
employing  a  telescoped,  small  FOV,  state-of-the-art 
interferometer  spectrometer  designs,  cold-filtered  detection 
bands  and  a  mosaic  focal  plane  of  background-limited  detectors. 
The  I FOV  and  telescope  sizes  were  suggested  by  the  Air  Force. 

The  selected  baselines  represent  different  levels  of  technology 
in  the  optics  and  focal  plane  designs, 

ju"rn"r„T  '  '  r--.  hi  n~ 

With  a  grating  equal  in  area  to  the  telescope  aperture  and  a  required 
resolving  power  of  50,000  (0.1  cm'1  at  5000  cm-*)  the  FTS  system  has  a 
throughput  advantage  of  30  for  the  assumed  IF0V  of  6.56xl0-i*  radians,1 


b. )  to  define  the  performance  (acnievable  noise-equivalent  radiance, 

NESR)  of  these  configurations  as  a  function  of  design/mission 
parameters  such  as  IFOV  and  available  measurement  time.  It  was 
decided  subsequently  to  identify  one  concept  as  the  baseline 
configuration,  to  evaluate  its  performance,  and  to  provide  nomo¬ 
graphs  that  allow  performance  scaling  with  respect  to  basic 
system  design  parameters  such  as  collector  diameter,  IFOV,  mea¬ 
surement  time,  f/number,  etc. 

c. )  to  define  a  spectral  data  processing  method  that  promises  near 

optimal  performance:  maximal  false  target  rejection  capability 
consistent  with  reasonable  detection  probability  and  small 
uncertainty  in  the  inferred  gaseous  amounts.  This  part  of  the 
study  analyzed  a  simple  technique  of  normalized  spectral  corre¬ 
lation,  and  yielded  prescriptions  for  evaluating  detection  and 
false  detection  probabilities  as  functions  of  system  NESR, 
amount  of  target  gas  present,  threshold  value  of  the  correlation 
coefficient,  etc.  The  computed  performance  predictions  indicate 
that  the  selected  method  is,  indeed,  very  close  to  optimum. 

d. )  to  determine  a  "minimum  detectable  quantity"  of  target  gas  for 

each  of  15  molecular  species  of  interest  to  the  Air  Force,  and 
the  uncertainty  in  the  deduced  amounts.  The  minimum  detectable 
quantity  may  be  used  to  determine  the  molecular  column  thickness 
that  will  yield  desired  levels  of  detection  and  false  detection 
probability  for  given  values  of  the  temperature  difference  be¬ 
tween  the  gas  and  the  underlying  earth  background.  This  part  of 
the  effort  also  determined  the  best  spectral  region  (absorption 
band)  for  the  detection  of  each  target  species.  These  initial 
estimates  of  minimum  detectable  quantity  and  uncertainty  ignore 
all  interference  effects  except  detector  noise,  and  assume  a 
particular  model  for  the  earth/atmosphere  scene. 

e. )  to  determine  how  normal  variations  in  the  scene  model  affect  the 

minimum  detectable  quantities. 

f. )  to  quantify  the  effects  on  minimum  detectable  quantity  of  other 

possible  interferences.  These  interferences  include  spatial 
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i  nhomogerte i  t  i es  in  the  atmosphere  or  background  (which  reduce 
the  effectiveness  of  the  background  subtraction  techniques), 
and  errors  in  the  predicted  atmospheric  spectral  transmittance 
and/or  target  gas  spectral  absorption  coefficient  (whose  product 
Tvagv  's  sPectral  correlation  reference). 

The  target  species  for  which  minimum  detectable  quantities  were  deter¬ 
mined  are  CH^,  Cl^l,  CO,  C02>  OF,  HBr ,  HC1,  HF,  HI,  HN03,  NH^  N20, 

N02  and  S02> 

The  minimum  detectable  quantity  (MDQ.)  of  target  gas  is  defined  as  the 
smal lest  value  of 

D  =  FuABv(Tg,Ts)  (1) 


that  can  be  detected  with  high  probability  (~95  percent)  and  low  false  detec¬ 
tion  probability  (~l  percent).  In  this  expression  F  is  the  fraction  of  the 
IFOV  occupied  by  the  target  gas  cloud,  u  is  the  column  thickness  of  the 
target  (molec/cm  ),  and  AB^  is  the  average  Planck  spectral  radiance  differ¬ 
ence  between  the  target  and  the  underlying  surface,  which  is  determined  by 
the  target  and  surface  temperatures,  Tg  and  Tg.  The  definition  of  D  and 
MDQ  derive  from  the  fact  that  the  observable  spectral  radiance  difference 
("contrast")  between  a  part  of  the  scene  containing  the  target  and  one  not 
containing  the  target  is,  ideally, 


AN 

v 


Dt  a 
v  gv 


(2) 


where  is  the  atmospheric  spectral  transmi ttance  between  the  target  and 
sensor,  and  agv  is  the  spectral  absorption  coefficient  of  the  target  gas . 

Detection  of  the  target  gas  cloud  in  an  actual  mission  would  consist 
of  first  computing  the  normalized  coefficient  of  spectral  correlation  p. 
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I 


defined  by 


P  = 


spectral  covariance  of  AN  and  r'o' 

_  v _ v  gv 

(standard  deviation  of  AN  )* (standard  deviation  of  r'a'  ) 

v  v  gv 


(3) 


where  *s  the  "reference  spectrum",  an  estimate  of  the  actual  spectrum 

Tvagv‘  lpl  exceeds  a  predetermined  threshold  R,  the  target  species 
would  be  presumed  present.  The  statistics  of  |p|  when  the  contrast 
spectral  radiance  ANv  contains  system  noise  (and/or  atmospheric/background 
interferences)  determines  the  probability  of  detection  P (T | T)  and  the  proba¬ 
bility  of  false  detection  P(TjO).  These  probabilities  vary  with  the  selec¬ 
ted  threshold  R,  the  number  of  resolved  spectral  elements  in  the  detection 
band,  the  variance  of  ANy  *=  Drvagv  relative  to  the  mean  square  spectrum 
noise  (NESR  ),  and  on  several  other  variances/covariances  that  characterize 
atmospher i c/background  i nterferences . 

The  amount  D  of  target  gas  present  is  inferred  from  a  least-squares 

fit  of  AN  to  the  reference  spectrum  x'a'  . 

v  v  gv 

The  basis  of  the  detection  and  quantification  methods,  and  the  proce¬ 
dures  or  codes  used  to  obtain  Tvaigv,  P(T[T),  P(T|0)  and  the  target  M0Q.‘s 
are  described  in  Section  2. 

The  baseline  system  configuration  used  in  computing  the  target  MDQ's 
is  described  in  Section  3  (which  refers  to  the  selected  concept  as  baseline 
Configuration  2).  This  system  configuration  uses  a  four-port  Michelson 
interferometer  with  two  dichroics,  which  allows  the  simultaneous  observation 
of  up  to  four  different  target  species.  The  fourth  port  (second  input  port) 
provides  for  illumination  of  the  back  side  of  the  interferometer  beamsplitter 
by  an  internal  blackbody  reference  or  by  a  different  portion  of  the  external 
scene  than  is  viewed  by  the  primary  input  port.  This  arrangement  effects  an 
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optical  subtraction  (gives  the  spectral  difference  of  the  displaced  scenes, 
or  the  scene  spectrum  minus  the  reference  blackbody  spectrum),  and  also  a 
potentially  large  reduction  in  the  interferogram  dynamic  range  (see  Vanasse 
et  alJ) .  The  basic  configuration  could  also  be  used  without  the  fourth 
port.  Then,  post-detection  subtraction  of  the  outputs  of  different  detector 
pairs  would  suppress  the  background  (perhaps  even  more  effectively  than  the 
four-port  optical  subtraction),  but  would  not  reduce  dynamic  range  in  the 
detector  outputs  ( interferograms) . 

Table  1.1  lists  the  gross  radiometric  parameters  that  enter  into  the 
computation  of  NESR  and  dynamic  range  for  the  selected  baseline  system.  The 


Table  1.1  Selected  Baseline  Configuration 


SYSTEM  FOUR-PORT  INTERFEROMETER  WITH  TWO  DJCHR01CS  AND  FOUR  SPECTRALLY 
OPTIMIZED  DETECTOR  ARRAYS 

Resolved  spectral  interval,  Av  ■  0.1  cm  '  (one-sided  interferogram) 
Spectral  range  2.3  -  13  pm 

Scan  time  (before  scaling)  ■  minimum  for  BLIP  performance 

OPTICS  Collector  diameter  D  ■  4  inches 
F#  *  3 

IFOV  «  0.656  x  0.656  mr  (“  10 m  footprint  at  50,000  foot  range) 
Overall  system  efficiency  £  *  0,195 
Emissivity  to  detector  c  ■  0.249 
Optics  transmission  t  ■  0.26) 


DETECTORS  For  3  to  4.8  pm  region, 

InSb,  n  -  0.5,  D*f*  -  5  x  1016 

For  7-0  to  12.7  pm  region,  assume 


HgCdTe,  n  ■  0.7,  D*f*  ■  8  x  lo'^ 

Element  size  x  »  9.5  mils  ■  0.24  mm 

Array  size  (four  arrays),  4  x  4  »  16  detectors  per  array 


selection  of  detector  materials,  quantum  efficiencies  n,  and  D*f*s  is  based 
on  Ref.  1.  Note  that  the  selection  of  HgCdTe  detectors  for  the  7.0-  12.7  um 
region  is  tentative;  extrinsic  Si:xx  detectors  (xx  »  Ga,  Bi  or  As)  are  also 
candidates  at  these  wavelengths.  Note,  also,  that  the  scan  time  is  speci¬ 
fied  as  the  minimum  corresponding  to  background- 1 imi ted  detector  performance; 
i.e.,  the  scan  time  for  each  detection  band  corresponds  to  the  frequency 
D  f  /BLIP  D  .  This  convention  is  used  merely  to  facilitate  scaling  of 
system  performance  according  to  a  pair  of  nomographs  (see  Section  1.4). 

Table  1.2  gives  the  emissivities  and  loss  factors  used  in  the  radio- 
metric  model  of  the  interferometer,  and  shows  how  the  effective  emissivity 
e*0.249  was  calculated.  Note  that  all  losses  are  represented  by  emissivi¬ 
ties,  and  that  obscurations  and  back-reflections  are  represented  as  elements 
of  unspecified  temperature.  An  example  is  the  secondary  mirror  and  its 
support  structure  (spider),  which  in  a  Cassegrain  telescope  system  is 
assumed  to  obscure  20  percent  of  the  primary.  The  spider  is  represented  as 
a  black  emitter  occupying  one  percent  of  the  clear  aperture.  The  primary 
and  secondary  are  assumed  to  have  a  combined  reflectance  of  95  percent,  or 
an  emissivity  of  five  percent.  The  dewar  window,  represented  by  two  ele¬ 
ments,  has  a  net  transmission,  after  reflection  and  absorption  losses,  of 
(1-0.08) (1-0.02)  a  0.902,  and  has  an  emissivity  of  0.02.  The  contribution 
of  each  300-degree  element  to  the  effective  system  emissivity  e  is  obtained 
by  multiplying  its  emissivity  by  the  product  of  transmi ttances  (1-e  values) 
of  the  elements  between  it  and  the  detector.  The  total  optics  transmission 
Is  simply  the  product  of  the  1-e  values  for  all  elements. 

The  computation  of  the  system  average  noise  equivalent  radiance  NESR 
for  each  detection  band,  and  the  corresponding  dynamic  range  in  the  inter- 
ferogram  is  described  in  Section  2.4.  These  two  performance  parameters 


Table  1.2  Radiometric  Mode)  of  Baseline  System 


Optical  Element 

Emisslvity 

c 

Temperature 

(*K) 

1 

Spider/secondary  obscuration 

0.20 

— 

2 

Spider  emission 

0.01 

300 

3 

Primary  and  secondary 

0.05 

300 

4 

Interferometer  (beamsplitter)  reflection 

0.50 

— 

5 

Beamsplitter  emission 

0.10 

300 

6 

Interferometer  mirrors 

0.05 

300 

7 

Dichroic 

0.10 

300 

8 

Dewar  window  loss 

0.08 

— 

9 

Dewar  window  emission 

0.02 

300 

10 

Cold  shield 

-- 

— 

11 

Cold  filter 

0.0 

-- 

12 

Detectors 

— 

(NA) 

(2)  Spider  0.01  x  0.3295  -  .0033 

(3)  Primary,  secondary  0.05  x  0.3469  ■  .0173 

(5)  Beamsplitter  0.10  x  0.7709  ■  .0771 

(6)  Interferometer  mirrors  0.05  x  0 . 8 l 1 A  ■  .0406 

(7)  Dichroic  0.10  x  0.9016  -  .0902 

(9)  Dewar  window  0.02  x  1 .0  ■  .02 

EFFECTIVE  EMISSIVITY  e  0.249 

TOTAL  OPTICS  TRANSMISSION  t  0.261 
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depend  on  the  detector  photon  flux  density,  in  addition  to  the  system 
radiometric  parameters  given  in  Tables  1.1  and  1.2  The  total  detector 
photon  flux  density  is  the  sum  of  the  system  internal  photon  flux  density 
to  the  detector  and  the  external  atmosphere-plus-background  photon  flux 
density  reaching  the  detector.  A  pre-release  version  of  the  AFGL  code 
FASC001  was  used  to  compute  the  atmosphere/background  spectral  radiances, 
which  were  converted  to  photon  spectral  radiances  and  then  integrated  to 
obtain  band  photon  radiances  (BPR's).  FASC0D1  was  also  used  to  obtain  the 
atmospheric  spectral  transmittance  and,  for  some  of  the  target  species, 
the  spectral  absorption  coefficient  a  ,  which  are  both  required  in  eval¬ 
uating  the  MDQ's  of  the  target  species  (see  Section  2.k). 

The  major  results  of  the  study  are  presented  in  a  standard  format, 
consisting  of  two  figures  and  one  table  for  each  candidate  detection  band 
(as  many  as  three  different  molecular  bends  or  segments  of  a  single  band 
were  considered  for  some  of  the  target  molecules).  Figure  1.1,  Table  1.3, 
and  Fig.  1.2  are  a  sample  of  the  results*,  they  represent  the  Vj  band  of 
SOj  (1090-1210  cm  *).  Results  in  this  standard  format  are  given  in  Vol .  II 
for  all  of  the  candidate  detection  bands.  All  of  the  tabular/plotted  values 
In  Vol.  II  are  based  on  the  Midlatitude  Summer  model  atmosphere  (see  Ref.  3 
and  Section  2.4.1),  and  assume  nadir  viewing  of  the  target.  Also,  they 
assume  no  atmospheric/background  interferences  (no  spatial  variations),  and 
that  the  reference  spectrum  x^a  is  known  exactly. 

The  uppermost  curve  of  Fig.  1,1  gives  the  spectral  radiance  between 
1090  and  1210  cm  ^  that  would  be  observed  by  a  nadir-viewing  sensor  viewing 
a  300K  earth  surface  through  a  Midlatitude  Summer  model  atmosphere. 

This  radiance  corresponds  to  the  case  of  no  SOj  in  the  detector  field-of- 
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view  (on  the  scale  of  the  figure,  the  spectral  radiance  added  by  S02  would 

be  barely  discernible,  since  the  minimum  detectable  quantity  is  very  small). 

The  middle  part  of  the  figure  shows  the  spectral  transmittance  of  the 

atmosphere.  FASCODI  was  used  to  obtain  the  radiance  and  transmittance 

spectra  and  also  to  degrade  them  to  a  spectral  resolution  of  0.1  cm'1 

(sin  v/v  Full  Width  at  Half  Maximum).  The  lowest  part  of  the  figure  shows 

the  "normalized  spectral  contrast"  of  the  $02  V|  band.  This  quantity  is 

a  x  normalized  to  its  maximum  value;  i.e.,  it  is  ideally  the  normalized 
gv  v  ' 

contrast  spectral  radiance  (Eq.  2)  and  the  normalized  correlation  reference 
spectrum. 

Table  1.3  (the  second  page  of  the  standard  output)  gives  various  param¬ 
eters  of  the  detection  band,  system  performance  parameters,  and  the  result¬ 
ing  minimum  detectable  quantity  corresponding  to  approximately  95  percent 
detection  probability  and  approximately  one  percent  false  detection 
probability.  Quantities  preceded  by  an  asterisk  would  change  if  any  of  the 
system  radiometric  parameters  were  changed.  The  spectral  resolution  should 
really  be  considered  a  constant  since  Av  ■  0.1  cm  *  is  close  to  an  optimal 
value.  However,  for  some  of  the  target  species  agv  had  to  be  determined 
from  measurements  corresponding  to  lower  resolution  (Av>0.1).  In  the 
standard- format  outputs,  the  MDQ  is  always  based  on  a  scan  time  equal  to 
"min  t^";  again,  this  convention  was  adopted  only  to  facilitate  scaling  of 
the  MDQ  results.  Table  1.3  gives  the  value  1  x  10^  for  the  MDQ  of  S02  (v^). 

It  should  be  noted  that  the  photon  fluxes  listed  in  the  summary  tables 
(e.g.,  Table  1.3)  do  not  Include  the  contribution  from  the  fourth  inter¬ 
ferometer  port.  The  fourth  port  approximately  doubles  the  listed  photon 

fluxes,  which  means  that  the  BUR  D*  should  be  divided  by  /2,  and  the 

c 

(NESR)qS  multiplied  by  The  effect  of  the  fourth  port  is  included, 


however,  in  the  listed  MDQ's  and  uncertainties  in  D.  In  other  words,  the 
tables  really  apply  to  the  case  where  background  subtraction  is  performed 
in  the  data  processing  (three-port  interferometer)  rather  than  optically 
(four-port  interferometer).  The  post-detection  subtraction  effectively  in¬ 
creases  the  NESR  by  /2;  hence,  both  techniques  yield  the  same  MDQ  and 
uncertainty.  The  fourth  port  would  be  used  primarily  to  reduce  dynamic 
range. 

Figure  1.2  displays,  for  the  SO2  band,  the  tradeoff  available  be¬ 
tween  detection  and  false  detection  probabilities  for  various  threshold 
levels  on  the  spectral  correlation  coefficient.  The  figure  shows  that  the 
threshold  value  R  corresponding  to  the  MDQ  given  in  Table  1.3  is  approxi¬ 
mately  0.07.  Note  that  the  figure  has  two  abscissa  scales.  One  is  used 
according  to  the  amount  D  of  target  gas  present.  The  other  gives  the 
corresponding  "correlation  s ignal-to-noise"  Z,  defined  as  o'D//2  NESR, 

where  a'  is  the  standard  deviation  of  the  reference  spectrum  x  a 

v  gv 

According  to  Eq.  (2),  Z  is  the  standard  deviation  of  the  observed  contrast 

spectral  radiance  AN^  divided  by  the  effective  rms  system  noise.  In  this 

example,  Z  is  approximately  0.12  at  the  MDQ  level.  This  means  that  the 

observed  AN  would  look  more  like  noise  than  the  characteristic  x  ct„  spec- 
v  v  gv  r 

tral  signature  of  the  SO2  band. 

It  is  worth  noting  that  Z,  at  the  MDQ  level,  is  always  between  approxi¬ 
mately  0.1  and  0.3  for  all  of  the  candidate  detection  bands  (see  Vol .  II). 
This  implies  that  it  is  generally  not  possible  to  determine  the  temperature 
of  the  target  gas  by  examination  of  the  distribution  of  rotational  line 
strengths  in  the  AN^  spectrum  (unless  D  >>  MDQ).  Thus,  it  is  not  possible 
to  determine  the  "minimum"  target  gas  column  thickness  u  independently  of 


AB(Tg,Ts)  solely  by  the  FTS  spectral  measurements.  However,  it  is  possible, 
for  many  of  the  detection  bands,  to  determine  the  background  temperature  T$ 
from  one  of  the  two  measurements  used  to  calculate  AN^,  provided  the  base¬ 
line  configuration  is  used  in  the  three-port  rather  than  (or  in  addition  to) 
the  four-port  mode. 

The  fact  that  the  target  gases  can  be  detected  and  quantified  using 
spectral  measurements  of  rather  low  s ignal -to-noi se  justifies  the  use  of  a 
spectral  correlation  process  that  involves  but  one  threshold  parameter.  In 
other  words,  the  use  of  a  multidimensional  decision  space  corresponding  to 
a  Bayes  test  (likelihood  ratios)  on  the  spectral  elements  would  not  be  ex¬ 
pected  to  yield  a  very  large  increase  in  detectability,  since  the  expected 
large  noise  tends  to  make  the  decision  region  a  hypersphere.  The  selected 
correlation  process  could,  perhaps,  be  made  more  nearly  optimum  by  the 
simple  procedure  of  excluding  the  minority  negative  or  positive  spectral 
contrast  radiances  AN^  (It  is  assumed  that  phase  correlation  would  be  used 
in  the  recovery  of  FTS  spectra). 

It  is  pointed  out  in  Section  5  that  spectral  elements  containing  strong 
atmospheric  absorptions  must  be  excluded  to  ensure  that  atmospheric/back¬ 
ground  spatial  variations  do  not  represent  overwhelming  interferences.  The 
fact  that  spectral  information  can  be  manipulated,  in  the  computer,  on  an 
element-by-element  basis  makes  the  proposed  technique  far  superior  to  the 
well-known  analog  method  of  gas-cell  correlation. 

The  performance  figures  given  in  the  standard  outputs  and  in  the  follow¬ 
ing  summary  tables  may  contain  small  errors  resulting  from  a  bug  in  the  pre¬ 
release  FASC001  that  was  not  discovered  until  most  of  the  computations  were 
completed.  The  FASC001  error  was  In  the  definition  of  function  F3,  which 


is  one  of  the  four  line  shape  functions  used  to  synthesize  the  Voigt  line 
shape  (cf.  Clough  and  Kneizys^).  This  error  can  have  a  significant  effect 
only  for  very  strong  lines  and/or  large  optical  depths,  and,  hence,  in  the 
present  application  is  important  only  when  the  target  species  is  overlapped 
by  very  strong  atmospheric  lines  (as  in  the  case  of  the  SC^  Vj  band,  shown 
in  Fig.  1.1).  Re-evaluation  of  some  of  the  computed  results  indicated  that 
the  NESR's,  MDQ's  and  computed  uncertainties  in  D  contain  errors  which  are 
generally  less  than  a  few  percent. 

1 . 1  Minimum  Detectable  Quantities 

The  following  tables  summarize  the  performance  capabilities  of  the 
selected  baseline  configuration  under  "standard"  and  "ideal"  conditions. 

The  tables  list  minimum  detectable  quantities,  and  also  the  (Corresponding 
NESRs.  The  MDQ's  scale  directly  as  the  system  NESRs,  and  the  required 
dynamic  range  is  inversely  proportional  to  NESR  under  conditions  of  BLIP 
performance.  The  complete  scaling  relations  for  NESR  and  dynamic  range  are 
provided  in  Section  1.4.  The  effects  of  non-standard  conditions  (climatic 
scene  variations)  are  summarized  in  Section  1.2.  The  effects  of  non-ideal 
conditions  (atmospheric/background  interferences  in  the  background- 
subtraction  process)  are  reported  in  Section  1.3. 

Table  1.4  lists  all  of  the  candidate  detection  bands  considered  in 
this  study,  the  computed  values  for  min  t^,  and  the  corresponding  (NESR) , 
MDQ's  and  uncertainties  in  the  minimum  detectable  quantities.  Note  that 
the  HC1  and  DF  bands  have  been  divided  arbitrarily  into  sub-bands  of  width 
100  or  200  cm  since  it  was  anticipated  that  detection  bands  of  width 
greater  than  200  cm  '  would  correspond  to  unrealistically  large  dynamic 


range  in  the  interferogram. 
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Table  1.5  gives  corresponding  system  performance  figures  scaled  to  an 

] 

_  I 

interferometer  scan  time  of  10  seconds.  In  the  BLIP  region,  system  NESR 

I/Z  ' 

scales  simply  as  [min  t^(sec)/10]  which  is  the  quantity  listed  as  ! 

"Scaling  Factor".  The  fifth  column  in  Table  1.5  is  the  minimum  detectable 
column  thickness  u,  in  molec/cm  ,  corresponding  to  a  temperature  difference 
of  1°C  between  the  target  gas  and  an  ideal  black  background.  Since  column 
thickness  (or  gas  density  or  mass  flow  rate)  is  assumed  to  be  of  greater 
interest  than  the  detectable  quantity  D,  we  used  the  former  minimum  quan¬ 
tity  to  select  the  "best"  detection  bands.  The  relationship  between 
MIN  u(AT=l°C)  and  MDQ.  is,  of  course,  strongly  dependent  on  the  spectral 
location  of  the  detection  band  (the  derivative  of  the  Planck  function  with 
respect  to  temperature).  The  underscored  MIN  u  values  in  column  5  indicate 
the  preferred  detection  bands. 

Table  1.6  lists  the  preferred  detection  bands  for  the  15  different 
target  species,  the  MIN  u's,  and  also  the  corresponding  minimum  detectable 
volume  concentrations  of  the  target  species  (in  parts  per  million  volume)  1] 

for  a  gas-cloud  thickness  of  10  meters  and  a  5°C  difference  in  iargwt-  ;j 

'  j 

background  temperature.  Note  that  the  COj  blue  spike  detection  band  is  a  j 

i : 

special  case,  corresponding  to  detection  on  the  basis  of  an  assumed  large  ' 

temperature  difference  rather  than  excess  column  thickness  (see  Section  2.4), 

Note,  also,  that  two  candidate  detection  bands  for  N0o  are  included  in  this 

/  l- 

! 

summary  table.  The  reason  is  that  the  preferred  NC^  band  requires  that  the  j 

system  include  a  detection  capability  extending  to  approximately  14.1  pm; 

i 

if  lower  performance  can  be  tolerated  for  NOj,  the  longwave  system  cutoff 
would  be  12.7  pm,  which  might  dictate  a  different  choice  of  detector  for 


the  longwave  detection  bands. 
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Table  1.5  Performance  Scaled  to  10-Second  Dwelt  Time  for  Candidate  Bands 


790-990 


Table  1.6  Seeled  Performance  for  "Best"  Bands 


Molecule 

Detection  Band  (cm  ') 

*MIN  u  (motec/cm2)  for  AT  -  1 *C 

**M1N  v  (ppmV) 

HF 

3240-3440 

2.6  *  1023 

2.1  x  106 

NCI 

2900-300C 

2.5  X  1018 

20 

N02(vj  +  Vj) 

2850-2935 

l.l  X  lo’9 

92 

OF 

2700-2900 

2.1  x  lO18 

17 

HBr 

2450-2650 

6.2  x  1 01 8 

51 

C02 

2375-2400 

(NA) 

*** 

0.0*47 

N20(v3) 

2160-2210 

2.6  x  lO18 

22 

CO(l-O) 

2130-2185 

1 .9  x  1C18 

16 

HI 

2100-2200 

6.5  x  1020 

5310 

ci W 

1295-1310 

2.0  x  lo’9 

163 

S02(v1 ) 

1090-1210 

2.0  x  lO18 

16 

NH3(v2) 

915-970 

1.8  x  lO17 

1.0 

HN03(2  vg) 

887.5-902.5 

2.0  x  JO17 

1 .6 

C2H4°2 

880-1000 

3.4  x  1017 

2.8 

ch3. 

79O-990 

3.9  x  tO18 

32 

N02(v2) 

710-795 

2.3  x  lo'8 

19 

SEE  following  page  for  footnotes. 
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Table  1.6  (Continued) 


* 

For  baseline  system,  except  dwell  time  t.  -  10  secs 

o 

* 

td  “  10  sacs.,  cloud  thickness  ■  10  m,  at  -  5*C,  T$fc  -  30OK 

MIN  v  »  WIN  u  (molec/cm2)  for  aT  -  i»c  6 

0.269  x  102  22i*£  *  L  (i000  cm)  x  6T(5*C)  * 

J  '  cur 

-  MIN  u  x  B. 17  x  ID'18 


For  td  -  10  secs.,  cloud  thickness  -  10  m,  and  aT  -  550*C  (T  -  85OK, 

Tsfc  "  J00K) :  9#S 

him  y  .  WIN  u  (molec/cm2)  for  flT  -  550*C  and  u  «  2400  cm'1  ,  6 
0.269  x  ,02  (jg)  x  ,000  X  '° 

(MIN  0  for  td  -  10  secs)/tB(850)  *B(300)  at  2k00  cm'1) 


0.269  x  I020  x  1000 


1  /  2 

(^)  2.5  x  IQn/(2.89  X  I0~4) 

0.269  x  I020  x  1000 


x  10° 


x  10° 


*  0.067  ppmV 
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The  minimum  quantities  given  for  N02  in  any  of  these  tables  may  be 
somewhat  pessimistic  because  of  the  model  used  for  the  vertical  distribution 
of  "normal"  atmospheric  N02»  The  atmospheric  vertical  column  thickness  for 
N02  was  approximately  1.7  times  greater  than  we  had  intended  it  to  be  as  the 
result  of  a  keypunch  error  (see  footnote  below  Table  2.1).  However,  the 
variability  and/or  uncertainty  in  atmospheric  N02  concentrations  is  probably 
on  the  order  of  a  factor  of  two.^ 

The  MDQ  and  MIN  u  for  HNO^  are  somewhat  pessimistic  because  the  detec¬ 
tion  band  includes  the  2  v(  transition,  but  not  the  nearby  vc  transition 

Zi  / 

(which  was  not  in  AFGL  trace  gas  line  parameters  compilation). 

Table  1.7  is  a  grtnd  summary  of  system  requirements  and  detection 
capabilities  for  the  15  molecules  of  interest.  The  table  segregates  the 
selected  shortwave  bands  (3“5  pm)  and  longwave  bands  (7.5-11*  urn) ,  and 
includes  the  dynamic  range  requirements  of  the  system.  The  listed  dynamic 
range  values  do  not  include  the  potentially  large  reduction  that  could  be 
effected  by  the  use  of  the  fourth  port  (second  input  port)  in  the  selected 
baseline  configuration. 


1 .2  Effects  of  Uniform  Scene  Variations 

The  MDQ's  will  depend,  obviously,  on  the  atmospheric  transmittance 
between  the  target  and  sensor,  and  also  on  the  brightness  of  the  scene, 
since  the  latter  influences  the  received  total  photon  flux  and,  hence,  the 
system  NESRs.  For  a  few  of  the  detection  bands  we  recalculated  the  MDQ's 
using  the  LOWTRAN  Midlatitude  Winter  model  in  place  of  the  Midlatitude 
Summer,  and  for  nadir  viewing  angles  equal  to  zero  and  l»5  degrees  (see 
Section  4).  Table  1.8  gives  the  lowest  and  highest  recalculated  MDQ's 


Table  1.7  Summary  of  Selected  Detection  Bends,  end  Performance  for  10-Second  Dwell  Time 


Table  1.8 

Variation  of  MDQ  Between  Midlatitude  Wi 
and  Midlatitude  Summer,  9-45° 

nter,  6»0, 

Molecule 

Detection  Band  (cm  ') 

MDQ  Relative 

(lowest) 

to  Standard 

(highest) 

HC1 

(2900-3000) 

0.73 

1.53 

N02(vj  +  v^) 

(2850-2935) 

0.75 

1.05 

DF 

(2700-2900) 

0.75 

1 .06 

HBr 

(2450-2650) 

0.74 

1.05 

CH4(va) 

(1295-1310) 

0.39 

2.80 

S02(v,) 

(1090-1210) 

0.53 

1 .26 

Midlatitude  Winter,  a  nadir  angle  of  zero  and  a  surface  temperature  of  273K, 
The  highest  correspond  to  Midlatitude  Summer,  nadir  angle  equal  to  45  degrees 
and  a  surface  temperature  of  300K. 

The  MDQ  variations  are  small  except  for  the  CH^  and  S02  bands.  The 
variation  for  these  two  bands  is  due  to  the  fact  that  they  contain  signifi¬ 
cant  HjO  absorptions,  and  that  approximately  five  times  more  water  vapor  is 
contained  in  the  inclined  summer  path  than  the  vertical  winter  path.  Over 
the  entire  range  of  latitudes,  seasons  and  zenith  angles  (0  to  45°),  the  H20 
path  column  thickness  could  vary  by  a  factqr  of  15. 

The  MDQ's  variations  given  in  Table  1,8  do  not  incorporate  any  kind  of 
geometrical  viewing  effects;  i.e,,  the  column  thickness  u  in  the  definition 
of  MDQ  is  measured  along  the  1 ine-of-s ight . 

1 .3  Effects  of  Spatial  Scene  Variations 

Section  5  deals  with  the  effects  of  scene  non-uniformity  (incomplete 


background  subtraction)  and  other  possible  types  of  interference  In  the 
detection  process.  A  major  interference  is  the  presence  of  differing 


atmospheric  HjO  amounts  in  the  two  parts  of  the  scene  that  are  differenced, 
since  this  can  result  in  very  strong  HgO  emission/absorption  lines  in  the 
contrast  radiance  spectrum.  In  fact,  we  show  in  Section  5  that  spectral 
resolution  elements  containing  any  strong  atmospheric  lines  (especially  H^O 
lines)  must  be  excluded  in  the  spectral  correlation  processing  in  order  to 
achieve  useful  MDQ  levels.  In  other  words,  the  effect  is  too  large  to  be 
tolerated  as  an  interference.  The  reduction  in  the  number  of  useable  spec¬ 
tral  elements  results,  of  course,  in  an  increase  in  the  MDQ. 

We  recalculated  the  MDQ's  to  reflect  this  modification  of  the  correla¬ 
tion  processing  only  for  the  selected  SOj,  HC 1  and  HBr  detection  bands.  The 
MDQ  for  SO2  increased  by  the  factor  1.5,  and  was  essentially  unchanged  for 
the  HC1  and  HBr  bands.  The  selected  SOj  and  HC1  bands  both  contain  strong 
HjO  lines,  but  there  is  less  coincidence  of  H2O  and  HCI  lines  than  there  is 
of  H20  and  S02  lines,  owing  to  the  very  different  nature  of  these  two  bands. 

The  effects  of  spatial  background  temperature  variations  are  also  con¬ 
sidered  in  Section  5,  but  only  for  the  HBr  detection  band.  It  is  shown  that 
a  difference  of  10°C  in  the  subtracted  backgrounds  would  lead  to  unacceptable 
performance,  unless  the  corresponding  Planck  spectral  radiance  difference 
in  the  contrast  spectrum  can  be  removed.  The  implication  is  that  the  con¬ 
trast  radiance  spectrum  must  be  "prewhitened"  before  it  is  correlated  with 
the  reference  spectrum. 

1 .4  The  Scaling  Nomographs 

The  performance  of  the  baseline  configuration  can  be  scaled  to  other 
scan  times  t^,  instantaneous  fields  of  view  a,  collector  diameters  0  and 
f/numbers  (F^) ,  using  the  summary  tables  of  Section  1.1  and  Figs.  1.3  and 
1.4.  Figure  1.3  is  a  plot  of  "Relative  NESR-otD"  versus  "Relative  t^";  the 
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RELATIVE 


Figure  1,4  Scaling  Nomograph  for  Detection  Band  Required  Dynamic  Range  (DR) 


RELATIVE 


point  (1,1)  represents  the  selected  baseline  configuration.  This  figure 
can  be  used  to  extrapolate  the  system  NESR.  The  detection  band  MDQ's  and 
uncertainties  in  D  scale  directly  as  the  NESR.  Figure  1.4  is  a  similar 
nomograph  for  scaling  the  system  dynamic  range  DR. 

The  region  of  background  limited  performance  is  a  single  line  in  each 
of  the  figures.  The  use  of  the  figures  outside  of  the  BLIP  region  (above 
the  lines  labeled  "BLIP")  depend  on  how  the  detector  D*F*  behaves  with 

^  ^  M  i  /a 

detector  area  A^.  For  InSb  detectors,  D  F  ~  while  for  S i : As  detec- 

^  ^  I  <2 

tors,  D  F  ~  .  The  nomographs  can  be  used  for  these  types  of  detectors 

^  -ff 

and  also  for  the  case  of  D  F  independent  of  detector  area.  There  appears 
to  be  no  simple  formula  relating  D  F  and  A^  for  HgCdTe.  Hence,  the  nomo¬ 
graphs  should  be  used  with  caution  for  the  detection  bands  that  use  HgCdTe 
detectors,  when  "Relative  td"  Is  close  to  unity  and/or  large  changes  are 
made  in  a,  D  or  F^.  The  following  examples  illustrate  the  use  of  the 
nomographs: 

EXAMPLE  A 

The  selected  dwell  time  is  made  equal  to  40  percent  of  min  t^  and  the 
product  Da  is  doubled  (e.g.,  a  might  be  doubled  while  D  is  held  fixed). 

iV  r~~~ 

The  system  has  a  detector  whose  0  F  scales  as  /A^;  hence,  we  use  the 
curve  labeled  "Relative  aD *  2  for  D*F*  ~  /A~J"  in  Fig.  1.3.  The  value  read 
from  the  abscissa  shows  that  NESR  x  aD  is  now  a  factor  of  two  greater.  But 
aD  was  doubled;  hence,  the  NESR  Is  unchanged.  System  performance  is  not 
in  the  BLIP  range;  also,  it  is  unaffected  by  changes  in  system  F^  for  this 
type  of  detector. 

Suppose  the  selected  dwell  time  is,  again,  40  percent  of  min  t^,  and 

£  A 

aD  is  unchanged.  Suppose  the  system  has  a  detector  with  constant  D  F 
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(independent  of  A^) .  If  we  change  the  F^  to  half  the  baseline  value,  the 
system  NESR  x  aD  doubles,  (we  use  the  curve  labeled  "Relative  F*  =  0.5  for 
OF  =  CONSTANT")  and,  hence,  the  NESR  doubles.  The  curve  we  use  for  this 
type  of  detector  is  not  dependent  on  changes  in  aD.  However,  if  aD  were 
doubled,  then  the  NESR  would  not  change  (since  NESR  x  aD  was  doubled). 

EXAMPLE  B 

Had  we  only  changed  the  dwell  time  to  0.40  x  min  t.  and  not  aD  for  the 

a 

system  whose  D  F  varies  as  the  NESR  would  increase  by  a  factor  of 

approximately  four.  Similarly,  had  we  only  changed  the  dwell  time  and  not 
Fy  or  aD  for  the  system  whose  D  F  is  constant,  the  NESR  would  also  increase 
by  a  factor  of  four. 

EXAMPLE  C 

Dwell  time  is  increased  to  four  times  min  td,  aD  is  halved,  and  the 

*  ?V  —  . 

detector  D  F  varies  as  vK^.  The  figure  shows  we  are  still  in  the  BLIP 
range  and  that  NESR  x  aD  is  half  its  original  value;  hence,  NESR  is 
unchanged.  Note  that  if  aD  were  not  changed,  we  would  use  the  same  curve, 
and  NESR  would  decrease  to  half  its  original  value.  Note,  also,  that  if 
dwell  time  were  increased  four-fold  and  aD  were  decreased  to  0.2  times  its 
original  value,  we  would  be  just  outside  the  BLIP  range;  this  would  also 

*  *  M 

happen  if  the  detector  D  F  were  constant  and  F  were  increased  to  five 
times  its  original  value. 

JL  JU 

The  nomograph  is  used  similarly  for  detectors  whose  D'F"  varies  as 

A^  (such  as  InSb),  except  that  the  choice  of  curve  in  the  non-BLIP 

#  2 

region  depends  on  the  relative  value  of  aD(F  )  .  The  dynamic  range  nomo¬ 
graph  is  used  in  the  same  fashion  as  the  NESR  nomograph. 
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2.1  Notation 


SECTION  2 


THEORETICAL  BASIS 


a.b.c 


Bv(T) 


B(Tg) 

B(T9) 


types  of  interferences 
area  of  collecting  aperture 
projected  area  of  gas  cloud 
area  of  detector 
footprint  area  of  sensor  1F0V 

Planck's  function  for  spectral  radiance  corresponding  to 
temperature  T 

Planck's  function  for  temperature  T9 

Planck's  function  for  temperature  T9 

Planck's  function  for  temperature  T9 

spectral  average  of  the  difference  6^(1^) -B^(Ts) 


AB  (T  ,T  )  spectral  average  of  the  difference  B  (T  ) -B  (T  ) 


v  g  s 


v  g  vs 


background-limited  performance 


Di  >  » 


BLIP  D' 


band  photon  radiance  (photons/cm  sr  s) 
speed  of  1 i ght 

=  hc/k:  second  radiation  constant 
diameter  of  collecting  aperture , 

ALSO  detectable  quantity  FuAB^(Tg,Ts)  of  target  gas 

detectable  quantities  of  other  (spectrally  interfering) 
target  gases 

band  detectivity  of  detector 
spectral  detectivity 
peak  spectral  detectivity 
background- 1 imi ted  detector  D 
dynamic  range  in  the  i nterferogram 
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f 


focal  length 

system  noise  bandwidth 


Af 

F  frequency  in  the  i nterferogram, 

ALSO  fraction  of  I FOV  occupied  by  target  gas 

F  maximum  frequency  for  BLIP  performance  (defined  by  detector 

.1  J- 

D  F‘  ) 
u 

F  f /number 

FTS  Fourier  transform  spectrometer 

G  G-factor  for  detector 

h  Planck's  constant 

H(R)  probability  distribution  function 

IFOV  instantaneous  field  of  view  (single  detector) 

J  photon  flux  density  on  detector 

JE  photon  flux  density  from  external  sources 

Jj.  photon  flux  density  from  internal  sources 

M  number  of  spectral  resolution  elements  in  detection  band 

MDQ  minimum  detectable  quantity  (minimum  D) 

MIN  v  minimum  detectable  volume  concentration  (ppmV)  for  t  .  = 

10  secs,  cloud  thickness  *  10  m,  Tg-Ts  =  5°C  and  Ts  =  300°K 

n  detector  noise  in  the  contrast  radiance  spectrum 

n  detector  noise  in  contrast  radiance  spectrum  and/or  spectral 

radiance  interferences  in  scene 

N^  spectral  radiance  of  scene 

AN  spectral  radiance  difference  between  two  1 FOV 1 s 

AN  spectral  radiance  difference  between  two  IFOV's 

v 

Nj  spectral  radiance  in  IF0V1 

N2  spectral  radiance  in  IF0V2 

atmospheric  spectral  emission  in  1F0V! 
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N*  or  NS 

N  S 

v 

NEP 

NER 

NESR 

NESR 


P 

P(T|T) 
P (T | 0) 
R 


Cd 

min  tj 
T 


T 

T 


g 

g 


T 

s 


T 


opt 


u 


X 


z  or  Z 

a 


a 


gv 


atmospheric  spectral  emission  in  1F0V2 

surface  radiance  in  IF0V1 

surface  radiance  in  IF0V2 

spectral  radiance  of  background  surface 

system  noise  equivalent  power 

system  noise  equivalent  radiance 

noise  equivalent  spectral  radiance 

average  noise  equivalent  spectral  radiance  over  detection 
band 

probability  density 
detection  probability 
false  detection  probability 

threshold  level  on  absolute  value  of  normalized  correlation 
coefficient  p 

dwell  time  (interferometer  scan  time) 
minimum  td  for  BLIP  performance 
temperature 

temperature  of  target  gas 
temperature  of  target  gas 
background  (surface)  temperature 
optics  temperature 

column  thickness  of  target  gas  along  the  1 ine-of-s ight 
=  c2v/(cT) 

unit  normal  random  variate 
detection  signal  to  noise 
angular  dimension  of  iFOV, 

ALSO  spectral  absorption  coefficient  of  a  target  gas 
spectral  absorption  coefficient  of  a  target  gas 
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spectral  absorption  coefficients  of  other  target  gases 
that  overlap  a 

estimate  of  a  or  a 

gv 

a  +  (D/Dj  )a  j  +  (D/D  2)a2  +  •  •  • 
or  (D/Dj)aj  +  (D/D2)a2 +  . . . 

effective  emissivity  of  the  optics  toward  the  detector 

quantum  efficiency  of  the  detector 

wavelength 

=  c/vp  longwave  cutoff  of  the  detector 
mean  of  AN 

v 

mean  of  i'a'  or  t  a 

v  gv 

mean  of  £/D 
wavenumber  (c/X) 

resolved  spectral  interval  or  "resolution" 
cAc 

detection  spectral  band  defined  by  cold  filter 

■  N2  ("no  gas")  -  N ^ :  difference  in  spectral  radiance  between 
two  IFOV's  when  no  target  gas  is  present  in  either  IFOV 

system  efficiency 

system  efficiency  due  to  non-optical  effects 

normalized  coefficient  of  spectral  correlation 

value  of  p  when  there  are  interferences  a  and/or  b  and/or 
c  but  no  system  noise 

standard  deviation  of  the  spectrum  to  =  to 

v  gv 

standard  deviation  of  the  spectrum  to  =  to 

v  gv 

standard  deviation  of  t'o' 
standard  deviation  of  £/D 
covariance  of  t6  and  x'a' 


°oc 

°IC 


T 


wd 


a 


covariance  of  tB  and  5/D 
covariance  of  -r'a'  and  5/D 
optics  transmission, 

ALSO  atmospheric  spectral  transmittance 

atmospheric  spectral  transmittance 

estimate  of  atmospheric  spectral  transmittance 

estimate  of  atmospheric  spectral  transmittance 

atmospheric  spectral  transmittance  in  IF0V1 

atmospheric  spectral  transmittance  in  IF0V2 

transmittance  computed  in  second  run  of  FASC0D1 

spectral  transmittance  of  target  gas  slab 

spectral  transmittance  of  target  gas  slab 

cosine  weighted  solid  angle  of  the  system  f-cone  at  the 
detector 

solid  angle  of  the  1F0V 


2.2  The  Observable  Quantity 

It  is  assumed  initially  that  all  detectors  (iFOV's)  observe  identical 
atmospheric  layers  and  background  surfaces,  i.e.,  identical  spectral 
radiances,  when  no  target  gases  are  present.  The  observed  spectral  radiance 
is  equal  to  the  atmospheric  spectral  emission  plus  a  background  contribution 

T  N  S  ,  (1) 

V  V 

where  tv  is  the  atmospheric  spectral  transmittance  and  N^S  is  the  intrinsic 
spectral  radiance  of  the  background  surface. 

A  slab  of  target  gas  of  spectral  transmittance  tj5  has  a  spectral  emis- 
sivity  of  I-tj*  and  an  intrinsic  spectral  radiance  equal  to  (I-t^)B^, 
where  =  B ^ (Tg )  is  the  Planck  spectral  radiance  corresponding  to  the 
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temperature  Tg  of  the  gas.  We  will  assume  the  target  slab,  lying  just  above 
the  surface,  is  very  thin;  i.e.,  that  the  atmospheric  transmittance  between 
the  top  of  the  slab  and  the  sensor  is  the  same  as  the  transmittance  of 
the  entire  atmosphere.  Then,  the  spectral  radiance  observed  in  an  IFOV  con¬ 
taining  the  target  gas  is  equal  to  the  atmospheric  spectral  emission  plus 


x  ( 1  -t  9)  B  9  +  t  t 

v  v  v  v  v 


9  N  S 


(2) 


The  target  gas  adds  to  the  observed  emission  via  the  first  term  and  attenu¬ 
ates  the  background  emission  via  the  second  term.  It  can  be  assumed  that 
the  observed  atmospheric  emission  is  not  affected  by  the  presence  of  a  very 
thin  target  slab. 

The  spectral  radiance  difference,  or  contrast,  between  an  IFOV  contain¬ 
ing  the  target  and  one  not  containing  the  target  is  the  difference  between 
expressions  (2)  and  (1): 


AN 

v 


T  (1-T9)  B9-  T  (  1  — T  9)  NS 
V  V  V  V  V  V 


(3) 


We  will  assume  the  target  gases  are  present  in  very  small  amounts.  That  is, 
the  slab  is  optically  thin,  which  implies 


Tv  ~  1  •  u  V  ’  (4) 

where  u  is  the  column  thickness  of  the  gas  along  the  1 ine-of-sight  and  o 

gv 

is  its  spectral  absorption  coefficient.  Then  the  contrast  spectral  radiance, 
from  Eqs.  (3)  and  (!*),  is 
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(5) 


AN 


T  a  u • (B  a 

v  gv  v 


NS) 

v 


We  will  make  the  additional  assumption  that  the  background  has  no  sharp 
spectral  structure  {emits  like  a  blackbody  of  temperature  T  ) ,  and  that  the 
detection  spectral  bands  are  narrow  {100  cm  '  or  less).  Then,  the  Planck 
radiance  difference  in  Eq.  (5)  can  be  approximated  by  its  spectral  average 
AB^  over  the  detection  band.  If  we  allow  for  the  occurrence  of  a  target 
that  does  not  fill  the  IFOV  (occupies  a  fraction  F  of  the  I FOV)  ,  the  ob¬ 
served  contrast  spectral  radiance  is 


AN  =  t  a  D  , 
v  v  gv 


D  -  Fu  4Bv(Tg,Ts)  . 


(6) 


D  is  defined  as  the  observable  or  detectable  quantity,  since  it  can  be 

determined  approximately  from  the  measurable  spectra  whose  difference  is 

AN  ,  and  a  prediction  of  t  a„  .  It  has  the  units  of  column  thickness  times 
v  v  gv 

2  2  -  ] 

spectral  radiance;  e.g.,  (molecules/cm  )(W/cm  srcm  ).  The  major  purpose 
of  this  study  was  to  determine  the  minimum  detectable  quantities  (minimum  D 
values  or  MDQ's)  of  specified  target  molecules  when  the  measured  AN^  includes 
a  noise  signal  representative  of  a  large-aperture,  state-of-the-art  FTS 
satellite  system.  In  principle,  the  high-resolution  spectral  measurements 
allow  an  independent  determination  of  the  gas  and  surface  temperatures; 
e.g.,  from  the  distribution  of  target  rotational  line  intensities  in  ANy  and 
from  the  envelopes  of  the  radiance  spectra  obtained  for  the  two  IFOV's. 
However,  by  using  the  correlation  method  described  in  Section  2.2,  it  is 
possible  to  detect  target  quantities  that  correspond  to  s ignal-to-noi se 
levels  much  less  than  unity.  In  other  words,  the  measured  contrast  spectrum 
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looks  very  much  like  noise  at  the  MDQ.  level  and,  hence,  it  cannot  be  used 

to  determine  line  intensities.  On  the  other  hand,  the  effective  surface 

temperature  can  be  obtained  easily  from  the  measured  spectra  for  one  or  more 

detection  bands.  The  implication  is  that  T  must  be  estimated  by  some  other 

9 

means  if  the  quantity  Fu  rather  than  D  is  required;  similarly,  an  indepen¬ 
dent  estimate  of  F  would  be  required  to  determine  u. 

The  fraction  F  is  the  projected  area  of  the  target  cloud  Ac  in  the 
IFOV  divided  by  the  known  IFOV  "footprint"  A  .  Thus, 

ArD  =  Acu  (7) 

■  (number  of  target  molecules  in  IFOV)  *  . 

In  other  words,  for  a  given  Planck  radiance  difference  (or  given  temperatures 
Tg’Ts)  the  sPectral  measurements  yield  the  number  of  target  molecules  in  the 
IFOV. 

Most  of  Section  2  deals  with  "the  ideal  case",  which  assumes 

a. )  the  atmospheric  spectral  transmittance  and  target  gas  spectral 

absorption  coefficient  a  can  be  predicted  exactly, 

b. )  the  absorptions  of  different  target  species  do  not  overlap  in 

the  observed  spectrum,  and 

c. )  the  atmospheres  and  backgrounds  in  the  two  IFOV's  are  identical. 
That  is,  system  noise  is  assumed  to  be  the  only  source  of  interference. 

Some  of  the  above  assumptions  are  relaxed  in  Section  5. 

2 .3  Method  of  Detection  and  Quantification 

The  measured  spectral  data  is  intended  to  provide  information  on  two 
I  eve  1 s : 

a.)  Does  the  data  imply  sufficiently  high  estimated  probability 
of  the  presence  of  a  target,  and  is  there  sufficiently  low 
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probability  that  noise  or  background  signals  will  be  falsely 
identified  as  a  target? 

b.)  If  the  above  conditions  are  true,  what  quantity  D  of  target 
is  present? 

Sections  2.3*1  and  2.3.2  address  the  first  question;  the  second  is  treated 
in  Section  2.3-3- 

2.3.1  Statistics  of  Detection/False  Detection 

One  must  first  prescribe  an  operation  on  the  spectral  data  that  will 
enable  recognition  of  the  target.  Generally,  this  operation  will  produce  a 
value(s)  to  be  compared  to  a  predetermined  threshold(s) .  If  the  threshold 
is  exceeded,  the  presence  of  a  target  is  declared.  The  fraction  of  the  time 
this  decision  is  correct  is  the  detection  probability  P ( T | T) .  The  fraction 
of  the  time  a  target  is  declared  when,  in  fact,  there  is  not  a  target  present 
is  the  false  detection  probability  P (T ( 0) . 

Generally,  by  decreasing  the  detection  threshold,  detection  probability 
can  be  gained  at  the  price  of  increased  false  detection  probability. 
Therefore,  in  selecting  the  threshold,  the  benefits  of  high  confidence  in  the 
target  event  must  be  balanced  against  the  risks  accruing  from  false 
identification.  A  benign  external  scene  and  low  system  noise  imply  that  high 
detection  probability  can  be  achieved  together  with  low  false  detection 
probabi 1 i ty . 

When  system  noise  and/or  external  interferences  are  included,  the  mea¬ 
sured  spectral  radiance  contrast  between  an  IFOV  containing  the  target  and 
one  not  containing  the  target  is 


0,  tv  and  a  were  defined  above;  is  a  spectrally  varying  quantity  that 
can  represent  detector  noise  and/or  the  effects  of  different  atmospheric  or 
background  properties  in  the  two  IFOV's.  A  near-optimal  process  for  target 
detection  is  based  on  the  spectral  correlation  between  aN^  and  a  prediction 
of  tvagv*  Specifically,  we  propose  to  use  (threshold  on)  the  magnitude  of 
the  normalized  correlation  coefficient 


where  denotes  a  spectral  sum  or  average  over  resolved  spectral  elements, 

u  is  the  mean  of  AN  and  u'  is  the  mean  of  x'a'  ,  which  is  a  theoretical 
v  v  gv 

estimate  of  the  actual  transmission-absorption  coefficient  product  T^otg^. 

Note  that  if  n^  =  0  in  Eq.  (8)  and  T^agv  =  Tvagv'  then  Negative  p 

corresponds  to  negative  D;  i.e.,  the  target  spectrum  is  observed  in  absorp¬ 
tion  because  the  target  gas  is  cooler  than  the  surface.  When  p  and  D  are 
positive  the  target  is  observed  in  emission.  Perfect  detection  corresponds 
to  |p|  =  l.  In  general,  |p|  is  less  than  unity  when  detector  noise  is 
present,  and/or  the  background/atmosphere  has  different  spectral  radiances 
in  the  two  IFOV’s. 

It  is  necessary  to  compute  the  probability  distribution  function  of  |p| 
in  order  to  determine  the  ability  of  the  system  to  detect  the  target  quantity 
D  and  reject  false  target  signals.  The  distribution  function,  H(R),  is  the 
probability  or  fraction  of  the  time  that  |p|  exceeds  the  value  R,  If  D  is 
not  zero,  this  is  the  same  as  the  detection  probability  corresponding  to  the 
particular  threshold  level  R  and  specified  value  of  D.  If  D  is  zero,  H(R) 
is  the  false  detection  probability. 


A  resonably  complete  determination  of  the  probability  distribution 
function  of  |p|  would  require  a  detailed  statistical  description  of  the 
spatial  and  spectral  properties  of  the  background/atmospheric  radiances, 


including  spatial  and  spectral  covariances.  This  information,  however,  is 

not  readily  available.  Our  only  recourse  in  determining  the  effect  on  p 

of  different  background/atmospheric  spectral  radiances  in  the  two  IFOV's  is 

to  evaluate  Eqs.  (8)  and  (9)  for  a  variety  of  different  severe  cases  or  a 

single  "worst"  case.  Computations  of  this  type  are  described  in  Section  5 

The  statistics  of  p  are,  of  course,  also  determined  by  the  detector  noise 

in  AN  . 
v 

We  will  treat  in  this  section  the  case  in  which  the  only  source  of 

error  in  the  measured  spectral  radiance  contrast  AN  is  detector  noise. 

v 

Then  the  n^'s  far  the  resolved  spectral  elements  can  be  regarded  as  independ¬ 
ent  random  normal  variates  with  zero  mean  and  standard  deviation  NESR, 
where  NESR  is  the  average  noise  equivalent  spectral  radiance  over  the  cold- 
filtered  spectral  band.  The  factor  results  from  the  fact  that  AN^  is  the 
difference  of  two  measurements  with  independent  noise  contributions.  We 

will  also  assume  j'a'  =  t  a  ;  i.e.,  that  the  atmospheric  transmittance  and 

v  gv  v  gv 

target  gas  absorption  coefficient  are  both  known  exactly. 

It  is  shown  in  Appendix  A  that  under  the  stated  conditions  (see 
Eqs.  (A-3)  and  (A- 13)), 
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where  o'  =  the  standard  deviation  of  the  reference  spectrum  t  a 

v  gv 

M  *  the  number  of  spectral  elements  in  the  summations 
Xj  ■  a  unit  normal  random  variate 

Z  -  o  'D/(vT  NESR)  -  the  ratio  of  rms  signal  to  rms  noise  in  the 
contrast  radiance  spectrum  (i.e.,  the  "detection  signal-to- 
noise") 


If  D  is  zero 
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Thus,  p  is  a  random  variate  dependent  on  M  random  noise  samples  n^  or  M 
values  selected  at  random  from  a  unit  normal  population. 

It  is  to  be  noted  that  according  to  the  second  form  of  Eq.  (10)  and 
Eq.  (11),  the  statistics  of  p  will  not  depend  on  the  detailed  shape  of  the 
spectrum  but  only  on  the  rms  variation  o'  of  the  spectrum  samples. 

This  is  illustrated  also  by  the  first  form  of  Eq.  (10)  when  it  is  noticed 
that  the  spectrum  samples  can  be  interchanged  in  any  manner  (since 

the  noise  samples  are  mutually  independent  and  uncorrelated  with  the  t  a 

v  gv 

samples)  without  affecting  the  value  of  p.  Thus,  for  example,  one  could 
reorder  the  samples  according  to  decreasing  amplitude  and  thereby 

convert  what  might  originally  be  a  very  structured  spectrum  into  a  rela¬ 
tively  smooth  one.  The  reordered  sample  set  would,  of  course,  have  the 


k2 


same  standard  deviation  as  the  original  set,  and,  therefore,  would  yield 
the  same  value  for  p. 


2.3.2  Evaluation  of  the  Detection  Probability  P(T|T)  and  False  Detection 


Probabi 1 i ty  P (T ( o) 

It  would  be  difficult  to  determine  analytically  the  probability  distri¬ 
bution  functions  of  the  variates  given  by  Eqs.  (10)  and  (11)  (the  P(T]T)  and 
P (T | 0) ,  respectively).  One  might  expect  that  closed  form  expressions  could 
be  obtained,  at  least  for  the  corresponding  probability  densities,  were  it 
not  for  the  fact  that  the  M  unit  normal  variates  in  the  numerators  are  the 
same  as  those  in  the  denominators;  i.e.,  the  numerator  cannot  be  regarded 
as  independent  of  the  denominator  unless  M  is  very  large.  We  elected  to 
evaluate  P(T|T)  and  P(T|0)  by  repeated  computations  of  Eqs.  (10)  and  (II) 
using  pseudo-random  numbers  generated  on  the  computer. 

Specifically,  we  first  assigned  a  value  to  M,  generated  M  random 
numbers  having  a  unit  normal  distribution,  and  computed  p(Eq.  10)  and 
p(Eq.  11).  In  the  case  of  p(Eq.  10)  it  is  necessary  to  perform  the  calcula¬ 
tion  for  various  values  of  Z.  This  process  was  repeated  5,000  times.  We 
then  determined  the  fractions  of  the  5.000  p(Eq.  10)  values  that  exceeded 
the  threshold  values  R  -  0.1,  0.15,  0.20 . 1.0.  These  fractions  repre¬ 

sent  P (T | T)  vs.  R  and  Z  for  given  M.  Similarly,  the  p  ( Eq .  11)  values 
were  used  to  obtain  P(Tj0)  vs.  R  for  given  M.  The  entire  procedure  was 
performed  for  several  values  of  M,  The  results  were  plotted  in  the  form 
P (T ) T)  vs.  P (T | 0) . 


*|f  the  numerator  and  denominator  of  Eq.  (11)  coutd  be  assumed  independent, 
p  (for  D  -  0)  would  obey  the  t-distrlbution,  which  approaches  a  normal 
distribution  for  large  M. 
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Figures  2.1,  2.2  and  2.3  represent  a  sample  of  the  computed  results. 

The  three  figures  correspond  to  signal-to-noise  values  Z  *  0.5,  1.0  and  1.5, 
respectively,  and  show  P (T j T )  vs.  P(TjO)  for  fixed  H  and/or  R,  Figure  2.2 
for  example,  shows  that  a  contrast  radiance  spectrum  with  20  resolved  spec¬ 
tral  elements  and  unit  signal-to-noise  Z,  yields  P(TjT)  =  93  percent  and 
P(T',0)  *  one  percent  if  the  threshold  on  |p|  is  0.55.  Alternately,  if 
93  percent  and  one  percent  are  the  required  detection  and  false  detection 
probabilities,  and  the  molecule  of  interest  is  observable  in  20  resolution 
elements,  we  would  be  required  to 

a. )  Set  the  threshold  on  |pj  at  0.55 

b. )  Achieve  Z  *>  1.0,  or  accept  a  detection  capability  0  corres¬ 

ponding  to  Z  *  o'D/(/2  NESR)  *  1.0. 

If  is  the  given  standard  deviation  of  the  20  spectral  elements  of  t  a  , 

v  gv 

and  NESR  denotes  the  achievable  noise  equivalent  spectral  radiance  of  a 
given  system  configuration,  then  the  minimum  detectable  quantity  MDQ.  is 

„dq  E  o  .  - 

min  a 

vT  NESR  ,. 

*  - ; -  (m  this  example) 

0 

The  three  figures  illustrate  the  expected  trends: 

a. )  Detection  performance  and  false  threat  rejection  capability 

both  increase  with  Z  and  rt. 

b. )  For  fixed  M  and  Z,  detection  and  false  detection  probabili¬ 

ties  both  decrease  with  increasing  threshold  R. 

What  is  surprising,  perhaps,  is  the  very  strong  dependence  of  system  per- 


'ornance  on  Z. 


Generally,  optimal  detection  performance  and  false  threat  rejection 
capability  will  be  realized  by  making  the  spectral  resolution  Av  comparable 
to  the  width  of  any  distinguishing  spectral  features  in  e.g.,  to 

the  line  halfwidth  if  the  spectrum  has  resolvable  line  structure  (then  Av 
should  be  approximately  0.1  cm  ').  That  is,  2  would  be  expected  to  have 
a  maximum  near  this  value  of  Av  because  of  the  way  o'  and  NESR  vary  with  Av. 
Note  that  as  Av  is  made  smaller,  o'  will  increase  asymptotically  to  a  con¬ 
stant  value,  while  NESR  «=  NER/av  will  continue  to  increase;  hence,  Z  will 
decrease.  If  Av  is  much  larger  than  the  width  of  spectral  features,  then 
generally  a'  will  decrease  more  rapidly  with  Av  than  NESR,  and  Z  will, 
again,  decrease.  Once  the  spectral  resolution  is  fixed  we  have  very  little 
control  over  the  values  of  o'  and  M  (if  there  are  two  or  more  IR  bands  of 
a  given  species,  we  would,  of  course,  choose  the  one  with  the  largest  o'  and/ 
or  M)  .  Thus,  in  using  Figs.  2. 1-2.3,  we  would  normally  regard  0'  and  M  as 
fixed.  If  the  system  performance  requirements  P(T|T)  and  P(T|0)  are  speci¬ 
fied,  and  the  NESR  has  been  calculated,  the  figures  would  be  used  to  deter¬ 
mine  the  minimum  detectable  quantity,  as  in  the  above  example.  If  the  MDQ 
so  determined  is  large  compared  to  D  values  expected  for  typical  targets, 
the  implication  is  that  the  system  must  be  redesigned;  e.g.,  given  a 
larger  collecting  aperture  to  achieve  a  smaller  NESR  and  proportionately 
sma 1 ler  MDQ. 

Volume  II  contains  plots  of  P(TjT)  vs.  Z  for  different  values  of 
P (T 1 0)  and  R,  for  each  of  the  candidate  detection  bands  (given  o'  and  M) , 
The  abscissa  also  have  D-scales  which  were  calculated  from  the  band  o's  and 
NESR's.  Plots  in  this  form  are  more  convenient  than  the  format  of  Figs,  2.1- 
2.3  for  scaling  the  MDQ  of  the  baseline  sensor  with  respect  to  NESR  (or 
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system  parameters  other  than  Av)  .  Like  Figs.  2. 1-2.3,  they  may  also  be 
used  to  trade  between  P(TjT)  and  P (T J 0) . 


2.3-3  Estimation  of  Qetectable  Quantity  D 

The  measured  spectral  radiance  contrast  is  ideally 

AN  =  0  t  a  (12) 

v  v  gv 


in  the  absence  of  noise.  An  optimal  method  of  estimating  D  from  the  spec¬ 
trum  aN^  is  to  find  the  value  D'  that  minimizes  the  mean  squared  spectral 
difference 
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(13) 


where  the  summation  is  over  the  M  resolved  spectral  elements  and  r'a'  is 

v  gv 

again  a  theoretical  estimate  of  x  a  .  Taking  the  derivative  of  c  with 

v  gv 

respect  to  our  estimate  D',  setting  the  result  to  zero,  and  solving  for  D' 
we  obtain 
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With  detector  noise  n  included,  the  spectral  radiance  contrast  is 

v 

AN  =  D  x  a  +  n  (15) 

v  v  g  v  v 

and  our  estimate  of  0  becomes 
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If  we  assume  t  and  a  are  accurately  known  O'a'  =  t  ot  },  then 
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This  estimate  is  merely  the  true  0  plus  a  weighted  sum  of  independent 
normal  variates  with  zero  mean  and  standard  deviation  Jl  NESR.  Using  the 
well-known  theorem  defining  the  statistics  of  a  sum  of  normal  variates 
(see  Appendix  A),  we  find  that  the  uncertainty  (standard  deviation)  of  our 
estimate  of  D  is 
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This  can  also  be  written 
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where  u'  and  o'  are  the  mean  and  standard  deviation  of  t  a  .  The  relative 
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uncertainty  is 


V 


771  [(’-)2*0)2] 


T72 


09) 


2 . ^  Computation  of  Normalized  Contrast  Spectrum  t  a 

■  —  ■  -  11  1  —  — ■  -  --r---  —  v  gv 

As  shown  above,  the  ability  to  detect  and  quantify  the  target  gases 

depends,  to  a  large  extent,  on  the  mean  and  variance  of  TvoigV>  which  will 

be  called  the  normalized  contrast  spectrum  (since  AN  equals  Dt  a  in  the 

v  v  gv 


50 


absence  of  noise  and  other  interferences),  or  the  "reference  spectrum", 
since  it  may  also  represent  a  predicted  spectrum  with  which  the  actual 

is  correlated.  The  target  gases  represented  by  the  spectral  absorption 
coefficients  a.^  are  of  two  different  types  in  this  study: 

a. )  Those  included  in  the  AFGL  atmospheric  and  trace  gas  line 

parameters  compilations. 

b. )  Other  molecules,  not  in  the  AFGL  compilations,  for  which 

the  Aerospace  Corporation  was  able  to  provide  ARC  with 

either  theoretical  calculations  of  a  or  spectral  measure- 

gv 

ments  of  this  quantity  (at  one-atmosphere  pressure  and 

temperatures  near  300°K) . 

The  "AFGL  molecules"  are  the  normal  constituents  of  the  atmosphere  and 
several  trace  molecules  present  at  all  levels  and/or  near  the  surface  in 
uncertain  amounts.  The  first  step  in  computing  tv<x  consisted  of  setting 
up  a  model  atmosphere  for  the  AFGL  molecules  (see  Section  2.^.1).  The 
vertical  distributions  of  temperature,  pressure  and  the  AFGL  molecules  with 
well-established  vertical  distributions  were  represented  by  the  AFGL 
(LOWTRAN)  Midlatitude  Summer  model  atmosphere.  The  distributions  of  some 
molecules  included  in  the  AFGL  line  parameters  compilations,  but  not  in  the 
LOWTRAN  model  atmospheres  (such  as  NO2  and  HNO^)  ,  were  based  on  information 
collected  by  the  Department  of  Transportation  (FAA)  during  the  HAPP  Program. 
AFGL  molecules  not  included  in  the  HAPP  data  or  LOWTRAN  model  (SO^  and  NH^) , 
and  the  "Aerospace  molecules",  were  assumed  to  exist  in  small  amounts  only 
near  ground  level,  as  localized  clouds. 

For  all  of  the  AFGL  target  molecules,  Tvagv  was  computed  on  the  basis 
of  two  runs  of  a  pre-release  version  of  AFGL's  FASC0D1  computer  code.  The 
first  run  obtained  the  spectral  transmittance  and  the  spectral  radiance 


of  a  vertical  path  through  the  adopted  atmospheric  model.  The  second  run 

obtained  the  spectra)  transmittance  (t  *)  under  the  same  conditions  except 

that  the  target  gas  density  in  the  bottom-most  layer  (0  to  10  meters)  was 

increased  by  a  nominal  ten  percent.  We  then  calculated  the  difference 

T  -t  *  and  divided  by  Au: 
v  v  ' 
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where  Au  denotes  the  excess  target  column  density  resulting  from  the  10  per¬ 
cent  increase.  The  term  t^(Au)  denotes  the  transmittance  of  the  excess 
column  density  of  the  taraet  gas  in  the  10-meter  slab,  which  is  always 
optically  thin.  Hence,  the  difference  tv-tv*  normalized  by  Au  reduces  to 

t  a  ,  which  is  independent  of  the  selected  Au  and  corresponds  to  an  a 
v  gv’  gv 

representing  ground-level  temperature  and  pressure. 

We  made  no  attempt  to  compute  the  target  gas  spectral  absorption 
coefficients  for  elevated  temperature  except  for  the  case  of  CO2  in  the 
i*.3  pm  "blue-spike"  region.  Increasing  the  temperature  would  generally 
cause  the  lines  to  become  narrower  (since  the  Lorentz  halfwidth  varies  as 
T  )  which  would  tend  to  increase  detection  capability.  It  would  also 
effect  a  redistribution  of  rotational  line  intensities.  However,  except 
for  the  case  of  molecular  hot  bands,  these  effects  would  be  expected  to  have 
an  unimportant  influence  on  the  MDQ.  That  is,  the  change  in  MDQ  would 
usually  be  small  compared  to  the  increase  in  the  quantity  D  available  for 
detection  caused  by  the  large  increase  in  aB^.  The  MDQ  need  not  be  known 
accurately  if  it  is  much  smaller  than  the  available  D. 
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In  the  case  of  CC^,  we  selected  the  detection  band  to  coincide  with 

the  blue-spike  region  of  the  4.3  pm  band.  Since  the  observable  blue-spike 

radiation  arises  from  a  hot  band  transition,  we  assumed  that  detection 

would  be  based  solely  on  a  large  temperature  increase  rather  than  an  excess 

column  thickness  u  and  a  small  AB^.  Accordingly,  the  two  FASC0D1  runs 

used  the  same  CC^  column  thickness  for  the  0-10  meter  layer,  but  different 

temperatures.  We  used  Eq.  (20)  to  obtain  ,  but  set  A u  equal  to  the 

total  CO,,  column  thickness  of  the  layer.  The  resultant  a  is  the  C0_ 

2  gv  2 

absorption  coefficient  at  the  elevated  temperature  (850°K)  used  in  the 

second  FASC001  run.  In  the  definition  of  MDQ  (Eq.  6),  u  then  represents 

the  column  thickness  of  the  at  elevated  temperature. 

In  the  case  of  molecules  for  which  Aerospace  provided  data  on  a  ,  we 

made  a  single  FASC0D1  run  to  obtain  t  ,  and  then  simply  computed  the  product 

t  a  of  the  merged  spectra.  The  t  a  spectrum  obtained  by  either  method 
v  gv  3  v  gv 

corresponds  to  a  spectral  resolution  (Full  Width  at  Half  Maximum)  of  approxi 
mately  0.10  cm  '  and  to  a  sinv/v  instrument  line  shape.  The  effect  of 
instrument  line  shape  was  introduced  in  the  FASC0D1  computations  via  its 
subroutine  SCANFN.  The  sinv/v  function  is  one  of  the  four  standard  options 
in  SCANFN. 

The  ARC  computer  code  CNTRST  used  to  obtain  t  a  from  the  FASC0D1/ 

v  gv 

Aerospace  spectra  also  obtained  the  standard  deviation  and  mean  of  T^a  . 

It  also  computed  the  net  photon  radiance  (by  integrating  the  photon  radiance 

spectrum  equivalent  to  the  FASC0D1  radiance  spectrum  obtained  with  t  ). 

The  photon  radiance  and  baseline  system  parameters  determine  the  system 

NESR  and  dynamic  range  (Section  2.5).  The  standard  deviation  and  mean  of 

t  a  ,  the  NESR,  and  the  P{T|T),  P(T|0)  nomographs  described  in 
v  gv 
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Section  2.3.2  determine  the  MDQ.  for  each  detection  band  and  the  uncertainty 
in  detected  MDQ. 

2.4.1  The  Atmospheric  Model 

The  atmospheric  model  used  in  most  of  the  computations  of  t  a  and 

v  gv 

L 

the  band  photon  radiance  is  the  LOWTRAN  Midlatitude  Summer  model.  This 
model,  and  the  five  others  in  the  LOWTRAN  codes,  give  vertical  distributions 
of  temperature,  pressure,  HjO,  CO2,  0^,  NjO,  CO,  CH^  and  O2 .  The  vertical 
distributions  used  for  NO,  NO2  and  HNO^  are  listed  in  Table  2.1  which  is 
based  on  measurement  data  reported  by  Luther^  (the  reported  data  shows  con¬ 
siderable  variations  between  measurements  from  different  sources;  the 
values  in  Table  2.1  represent  approximate  mean  distributions).  The  molecules 

S02  and  NH^  were  assumed  to  be  present  only  in  the  bottom  10-meter  layer  in 
1 3  2 

the  amount  1  x  >0  molec/cm  .  The  AFGL  code  ATMPTH  was  used  to  integrate 
the  LOWTRAN  Midlatitude  Summer  model  (to  find  mean  temperatures  and  pressures 
and  the  column  densities  of  the  species  1^0,  C02,  0^,  N20,  CO,  CH^  and  02  for 
the  15  layers  indicated  in  Table  2.1). 

2 . 5  Computation  of  the  System  NCSR  and  Required  Dynamic  Range 

This  section  gives  the  equations  used  to  obtain  the  background- 1 i mi  ted 
specific  detectivities  (BLIP  D*s)  of  the  system  detectors  and  the  corres¬ 
ponding  system  NESR's  and  dynamic  range.  It  is  emphasized  that  these  results 
assume  the  detectors  are  operating  at  frequencies  where  BLIP  performance  is 
achievable.  The  variation  of  D*  and  NESR  with  frequency  outside  of  the  BLIP 
range  is  discussed  in  Section  2,6. 

For  an  ideal  photoconduct i ve  BLIP  detector,  the  spectral  detectivity 
is  linear  with  wavelength  out  to  its  longwave  cutoff  Xc: 


5 


Table  2.1  Vertical  Distributions  used  for  NO,  NOj  and  HNO^ 

LAYER  COLUMN  THICKNESSES  (molec/cm2) 

Layer  Boundaries 


(km) 

NO 

N02 

HN03 

0-0.01 

3.9 

X 

1012 

1  .0  x 

io13 

1  2 

1 .0  x  10 

0.01-4.0 

1 .24 

X 

io'5 

3.1  X 

,o'5 

4.4  x  IO14 

4-7 

6.6 

X 

io'4 

1  . 15  X 

io15 

4.2  x  IO14 

7-H 

6.4 

X 

,o14 

9.0  x 

io’4 

7.2  x  IO14 

11-14 

3.0 

X 

io14 

5.1  x 

io14 

7.2  x  IO14 

14-17 

2. 1 

X 

io14 

4.8  x 

io'4 

9.6  x  IO14 

17-20 

1.5 

X 

,o14 

5.15  x 

io'4 

1 .2  x  IO15 

20-25 

3.0 

X 

io14 

9.9  x 

io14 

1 .85  x  IO15 

25-30 

4.25 

X 

io'4 

1  .0  x 

io'5 

1  k 

7.0  X  IO1 

30-35 

3-75 

X 

io14 

7.5  x 

1  5* 

10  * 

14 

1 . 15  x  10 

35-40 

2.1 

X 

io'4 

2.5  X 

io14 

2.2  x  1 O1 3 

40-45 

1.1 

X 

io’4 

1  .0  x 

io14 

0 

45-54 

0 

7.0  x 

io13 

0 

54-75 

0 

0 

0 

75-® 

0 

0 

0 

^  I  *4 

This  value  was  used  in  the  computations;  it  should  have  been  7-5  x  10 
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1 


where  the  peak  spectral  detectivity  D"  is  given  by 


!/ 2. 

n  X, 


1/2 


2hcJ 


1/2 


2hv  (Jp  +  J  ) 

c  t  i 
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(21) 


(22) 


where  n  is  the  detector  quantum  efficiency,  h  is  Planck's  constant,  c  is 
the  speed  of  light,  and  J  is  the  total  photon  flux  density  incident  on  the 
detector.  The  contribution  to  J  from  sources  external  to  the  sensor  is 
obtained  from 


flAt  fV 2  , 

aT  J  i^dv 

d  Vj 


TIT 


(2F#)2  v, 


v,  N 

/  2^dv 


hv 


(23) 


where 


2 

n  =  solid  angle  IFOV  of  the  system  =  a  ,  the  angular  IFOV. 

A  =  collecting  aperture  area 
t  =  optics  transmittance 
A^j  *  detector  area 
F^  *=  system  f/number  at  detector 

■  external  radiance  spectrum  (computed  using  FASC0D1) 
vj,v2  *  spectral  band  limits  defined  by  the  cold  filter  (v2  <_  vc) . 


The  second  form  of  Eq.  (23)  is  obtained  from  the  first  by  substituting 
2  o  a 

A  *=  nO  A,  ■  f  fi  and  Fff  *  f/D,  where  f  is  the  system  effective  focal 
length  and  D  is  the  effective  aperture  diameter.  (Throughout  this  report, 
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D,  and  or  a  =  Jv.  are  regarded  as  the  fundamental  geometric  parameters 
of  the  system.) 

The  contribution  to  J  from  internal  sources  (optics  emission  within 
the  coiu  shielded  f/cone)  is  given  by 


u  2 

where  =  ir/[I  +  (2Fff)  ]  is  the  cosine-weighted  solid  angle  of  the  system 
f-cone  at  the  detector,  e  is  the  effective  emissivity  of  the  optics  towards 
the  detector,  B^fT)  is  the  Planck  radiance  spectrum,  T  =  300°K,  x  =  c.,v/cT 
and  c2  *  hc/k,  the  second  radiation  constant.  The  series  in  Eq .  (2*0  (of 
which  only  the  first  five  terms  are  used  in  this  study)  is  the  infinite 
series  expansion  of  a  finite  integral  of  the  Planck  spectral  photon  emittance 
of  a  blackbody.^  The  series  converges  most  rapidly  for  large  x  (small  AT), 
and  five  terms  gives  sufficient  accuracy  for  wavelengths  as  long  as  30  urn 
when  T  =  300°K. 

The  system  performance  parameters,  such  as  noise  equivalent  radiance, 
depend  on  an  average  or  band  D*'  which  is  related  to  the  peak  spectral  detec- 

J. 

tivity  D"  through  the  so-called  G  factor: 
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where 


where  N^  is,  again,  the  external  radiance  spectrum.  This  equation  is  anal¬ 
ogous  to  expressions  which  relate  blackbody  and  monochromatic  detectivities 
Note  that  the  only  system  parameters  on  which  J  *  depends  are 

system  f/number,  optics  t  and  e,  and,  of  course,  the  two  spectral  band 

cutoffs.  Consequently,  and  D  (vj  to  v^)  depend  only  on  these  system 

C  H 

parameters  and  on  the  detector  quantum  efficiency;  i.e.,  on  F  ,  e,  t,  Vj, 

\>2  and  n. 

The  system  noise  equivalent  radiance  at  the  entrance  aperture  can  be 
written  in  the  alternate  forms 


,rJ/2 

I 

n  0* (v^  to  v^) 

F#Uf)'^2  _1_ 

X-j-OQ^  0  (vj  to  V2) 


(27) 


where  is f  =  n/C^t^)  is  the  noise  bandwidth,  t^  is  the  dwell  time  (interfer* 
ometer  scan  time),  and  X  Is  the  net  system  efficiency,  equal  to  r£e,  where 
5e  is  the  efficiency  due  to  non-optica!  effects  (e.g.,  the  signal  modulation 
efficiency) . 
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The  average  noise  equivalent  spectral  radiance  of  a  multiplexing 
spectrometer  (e.g.,  an  interferometer)  is  equal  to 


NESR  =  NER/Av  , 


where  Av  is  the  resolved  spectral  interval.  The  actual  NESR  is  given  by 


4-  NESR 


u,  +v2 


2y  NER 

v  |  +  \>2  Av 


F*(Af) 


2v  1 


(n,1+v,2Uv  D*(v,  to  v2) 


where  v  in  the  first  form  is  (vj+v2)/2,  Equation  (29)  results  from  the 
fact  that  the  spectral  NER  is  inversely  proportional  to  the  spectral  D" 
which  is  inversely  proportional  to  the  wavenumber  v;  i.e,,  NESR  is  directly 
proportional  to  v  and  its  average  over  (vj  ,\>2)  must  be  equal  to  NESR. 

The  required  dynamic  range  in  the  i nterferogram  defined  in  terms  of 
NER  and  other  quantities  referred  to  the  aperture  is 


Dynamic  Range  (DR)  = 


,  T  v_ 

—■  J  2  N  dv  +  — 
1ER  1  v  t 

Lvi 


/  B  (T  ) dv 
1  v  opt 

V1 


where  the  second  term  in  brackets  is  the  "signal"  contribution  of  optics 
emission  after  extrapolation  to  the  aperture.  In  terms  of  photon  flux 
density  at  the  detector,  this  can  be  written  (using  Eqs.  25,  26  and  22-2A) 


FW/2n’/2E  J,/2 


DR  -  t 
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Note  that  the  factor  J 


in  Eq.  (31)  is  approximately  proportional 


1/2 

to  1/Fff  (see  Eqs.  23  and  2*0.  The  approximation  assumes  that  the  denominator 
factor  1  +  (2F*)^  can  be  replaced  by  (2F^)^  or  that  J^>>JT.  This  means 
trac  the  dynamic  range  is  very  nearly  independent  of  f/number  for  a  cold- 
shielded  system  with  fixed  ft  and  D. 

The  equations  given  above  were  programmed  for  the  computer  to  obtain, 
for  the  various  detection  bands,  the  peak  spectral  detectivity  within  the 
passband  (D^  ),  the  band  detectivity  D  (v^  to  V2) .  the  average  noise 
equivalent  spectral  radiance  NESR,  and  the  system  dynamic  range,  The 
computed  values  apply  to  the  selected  baseline  system  and  are  reported  in 
Sect i on  1  and  Vol  .  II. 


2 . 6  Basis  of  the  Scaling  Nomographs 

The  nomographs  presented  in  Section  l A  may  be  used  to  scale  the  NESR 
(MOQ's)  and  requiicu  dynamic  range  of  the  selected  baseline  system  with 
respect  to  its  basic  radiometric  parameters,  such  as  collector  diameter, 
IFOV,  and  scan  time.  This  section  defines  the  construction  of  the  two 
scaling  nomographs. 

According  to  Eqs.  (27)  and  (28)  of  the  previous  section,  NESR  is  pro¬ 
portional  to 


NESR 


fW/2 

d//2  ; 

F*  I 

T7T^ 


Da  t 


(32) 


if  the  system  resolution  (Av)  and  efficiency  F,  are  fixed.  But  to  a  very 
good  approximation 
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'■N/ 


(33) 


F* 


for  a  cold-shielded  detector,  if  the  detector  is  operating  in  the  BLIP 

domain  and  c,  t,  n,  and  T  .  are  fixed  (see  Appendix  B) .  Since  D"  ~  D'  , 

o  p  t  I  c  s  ^ 

c 

we  have  from  Eqs .  (32)  and  (33) 


NESR 


(aD)  t 


1/2 


in  the  BLIP  domain. 

The  sketch  below  shows  schematically  the  variation  of  system  D  with 

A 

C 

electrical  frequency  F.  For  sufficiently  low  F  the  system  will  achieve 
BLIP  performance,  while,  for  higher  frequencies  (F  >  F")  ,  performance  will 
be  limited  by  detector-preamplifier  noise. 


The  limit  on  due  to  detector-amplifier  noise  will  vary  as 
c 


(35) 
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1 


‘1 

1 

1 


2 

for  some  types  of  detector;  e.g.,  an  arsenic-doped  silicon  detector.  The 
0  depends  on  detector  area  because  this  type  of  detector  is  more  properly 
characterized  by  an  area- independent  NEP  (noise-equivalent  power)  rather 
than  a  D  . 

The  highest  electrical  frequency  at  which  the  detector  must  operate 
is  the  nominal  number  of  interferogram  samples  V2/AV  divided  by  the  scan 
time  t^: 


F 


v2 

tdAv 


(36) 


Substituting  (36)  into  (35)  gives 


F^ (aD) td  , 


and  substituting  (37) 


into  (32),  with  O'  ~  O’,  gives 


NESR 


(37) 


(38a) 


for  the  high-frequency  (non-SLIP)  domain.  The  equivalent  relation  for 
detectors  with  area- i ndependent  D*f’’  (i.e.,  D*  ~  1/r)  is 

A 


NESR 


(aO)t 
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(38b) 


*c 

For  an  InSb  detector,  whose  OF  is  proportional  to  A 


-1/2 


,  the  equivalent 


NESR 


(38c) 


Equations  (3*0  and  (38a)-(38c)  are  the  desired  scaling  relations  for 
NESR  for  the  BLIP  and  high-frequency  domains,  respectively. 

A  convenient  nomograph  based  on  the  above  scaling  relations  is  shown 
schematically  in  Fig.  2.4.  The  quantity  plotted  is  the  relative  NESR  x  aD 
vs.  the  relative  dwell  time  td .  The  point  (1,1)  denotes  the  selected  base¬ 
line  system,  whose  dwell  time  is  defined  to  be  mint^,  the  lowest  value  for 
which  BLIP  performance  can  be  achieved.  According  to  Eqs.  (3*0  and  (38), 

NESR  x  «D  varies  td  for  td  1  s  larger  than  mint^  (BLIP  domain),  and  as 
-3/2 

tj  for  relative  td  1  s  less  than  unity.  Hence,  the  log-log  plot  has  two 
lines  of  slopes  -1/2  and  -3/2  through  the  reference  point  (1,1);  these 
curves  are  applicable  provided  only  td  is  varied.  The  other  lines  of  slope 
-3/2  in  the  non-BLIP  region  describe  the  NESR  x  aD  variations  with  aD  and/or 
for  the  three  types  of  detectors.  Note  that  only  one  set  of  labels  for 
these  curves  is  applicable  in  a  specific  case,  depending  on  which  of  the 
three  types  of  detector  is  employed. 

From  Eqs.  (31)  and  (22)  of  Section  2.5,  the  dynamic  range  (DR)  in  the 
interferogram  varies  as 


DR 


r'oDJ'/2  F#«'>'d'/2 

“w?77'  ~  o* 


£  u  9 

But  in  the  BLIP  domain,  D  ~  F  to  a  very  good  approximation,  and 


DR  ~  (aD) t 


1/2 


(39) 


(40) 
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*>-■ 

f 


x  1 


In  the  region  of  detector -amp li f i er  limited  performance  D"  varies  as 


Eq.  (37)  if  DV'  -  A, 


1/2 


Combining  (37)  and  (39),  we  obtain 


~  7772  ■ 

t  A 


(Ala) 


(again,  Av  is  assumed  fixed).  The  equivalent  relation  for  detectors  with 


k  _k 


area- independent  OF  is 


OR 


(<»D)Fff 

TT7F 


(Alb) 


while  for  detectors  wi th  D  F  ~  A^  ,  it  is 


DR 


(aD)2(F*)2 


(Ale) 


A  nomograph  based  on  the  dynamic  range  scaling  relations  is  shown 
schematically  in  Fig.  2.5.  This  nomograph  is  used  in  the  same  manner  as 
the  NESR  nomograph  of  Fig.  2. A.  More  detailed  versions  of  the  nomographs 
and  examples  of  their  use  are  given  in  Section  T ,  A . 
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SECTION  3 


THE  BASELINE  SYSTEMS 

Early  in  the  study,  we  developed  four  FTS  system  concepts  representing 
levels  of  technology  ranging  from  "completely  proven"  to  beyond  the  state- 

X 

of-the-art.  The  original  intention  was  to  determine  the  detection  capa¬ 
bilities  of  four  corresponding  configurations  for  a  variety  of  target, 
atmospheric  and  background  conditions,  as  a  function  of  target  range  and 
available  measurement  time. 

It  was  realized  subsequently  that  the  performance  of  the  most  advanced 
concept  could  not  be  estimated  unless  efforts  beyond  the  scope  of  the  study 
were  undertaken  to  determine  the  characteristics  of  the  large-mosaic  detectors 
and  some  of  the  optical  elements.  Also,  the  consensus  of  those  present  at 
an  early  Air  Force  program  review  meeting  was  that  one  of  the  proposed  con¬ 
cepts  was  basically  not  viable,  but  had  features  that  could  be  usefully 
employed  in  one  of  the  less  advanced  concepts.  Thus,  the  number  of  baseline 
configurations  was  reduced  to  three,  and  the  decision  was  made  to  evaluate 
only  the  second,  and  not  discuss  the  third  further.  It  was  agreed  also 
that  nomographs  would  be  developed  to  allow  the  performance  of  the  selected 
baseline  configuration  to  be  scaled  with  respect  to  basic  radiometric  param¬ 
eters  such  as  collector  diameter,  detector  I FOV  and  measurement  time. 

Tables  3-1  and  3.2  define  the  general  characteristics  of  the  base¬ 
line  FTS  configurations.  Configuration  1  consists  of  a  conventional  Michel- 
son  interferometer  cube  and  drive  mechanism  with  a  four  by  four  mosaic  of 
detectors,  a  series  of  cold  filters  that  can  be  introduced  sequentially  (or 

A 

The  system  concepts  were  selected  with  the  help  of  Bartlett  Systems,  Inc. 


Table  3.1  Baseline  Configuration  i 


SYSTEM  CONVENTIONAL  SINGLE-PORT  MICHELSON  INTERFEROMETER  WITH  MOSAIC 

FOCAL  PLANE 

Resolved  spectra!  interval,  Av  =  0.10  err,  '  (one-sided 
i nterferogram) 

Spectral  range,  3-3-12.6  pm  (one  beamsplitter) 

Number  of  cold  filtered  spectral  bands  —  20 
Cold  filter  spectral  width  ~  25  cm-'  to~200  cm-' 

Scan  time,  TBD  for  each  species;  nominally  10  secs. 

Scans  per  target  (number  of  bands  per  target  facility  or 
co-adds  per  single  band),  10 


OPTICS  Collector  diameter,  0  »  4  inches 

F#  -  3 

I F0V  =  0.656  mr  (2  10  m  footprint  at  50,000  foot  range) 
Overall  system  efficiency,  X  —  0.22 
Emissivity  to  detector,  c  cs  0.17 
Optical  transmission,  t  s:  0.29 

Cooling:  cold  shield,  spectral  filter  and  f-cone  cooled 
sufficiently  to  contribute  negligible  background. 
Telescope  uncooled  or  cooled  passively. 


DETECTORS  Material:  Si:Ga,  Si :Bi  or  HgCdTe 

Element  size,  x  =*  9.5  mils  *  0.24  mm 
Array  size,  4  x  4  =  16  detectors 


signal 

PROCESSING  Subtraction  of  pairs  of  IFOV's 
Gain  switching,  6  bits, 

A/D  conversion,  15  bits  (32,000  resolution  levels) 


Table  3.2  Baseline  Configuration  2 


SYSTEM 


OPTICS 


OETECTORS 


SIGNAL 

PROCESSING 


ADVANTAGES 


FOUR-PORT  INTERFEROMETER  WITH  TWO  DICHROICS  AND  FOUR 
SPECTRALLY  OPTIMIZED  DETECTOR  ARRAYS 

Resolved  spectral  interval,  Av  =  0.10  cm  '  (one-sided 
interferogram) 

Spectral  range,  3. 3" I  2. 6  ym  (one  beamsplitter) 

Number  of  cold  filtered  spectral  bands, o:  20 
Cold  filter  spectral  width  ~  25  cm"'  to~200  cm”' 

Scan  time,  t^  =  minimum  value  for  BLIP  performance. 

Scans  per  target  (number  of  bands  per  target  facility  or 
co-adds  per  single  band),  40 


Collector  diameter,  D  *  4  inches 
F  #  =  3 

IFOV  =  0.656  mr  (o-|0  m  footprint  at  50,000  foot  range) 
Overall  system  efficiency,  £  0.20 

Emissivity  to  detector,  e  =*  0.25 
Optical  transmission,  t  ^  0.26 

Cooling:  cold  shield,  spectral  filter  and  f~cone  cooled 
sufficiently  to  contribute  negligible  background. 
Telescope  uncooled  or  cooled  passively. 

Optical  layout:  see  Fig.  3.1 


For  3 -2-3 -7  ym  region,  InSb 
For  3- 7-**. 5  ym  region) 

For  4. 5-4.8  ym  region/ 

For  7.5-13.0  ym  region,  HgCdTe  or  Si:xx  (xx  =  Ga,  Bi  or  As) 
Element  size,  x  =  9-5  mils  =  0.24  mm 

Array  size  (four  arrays),  4  x  4  =  16  detectors  per  array 


Background  reduction  effected  by  second  input  port  (2nd  F0V 
or  internal  reference).  May  not  require  gain  switching. 


Dichroic  separation  permits  simultaneous  examination  of  four 
spectral  regions  (simultaneous  detection  of  more  species, 
and  greater  detectability  of  species).  Significant  dynami 
range  reduction  due  to  second  input  port. 
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possibly  a  cold  CVF) .  The  system  operates  at  a  spectral  resolution  of 
0.1  cm  *,  and  has  an  IFOV  (each  detector)  of  0.656  mr.  The  scan  time  is 
considered  a  variable  that  can  be  adjusted  according  to  mission  constraints 
and/or  the  number  of  different  target  species  that  must  be  detected  during 
the  available  measurement  time.  It  was  known  (from  Ref.  2)  that  the  maxi¬ 
mum  dynamic  range  in  the  i n terferog ram  would  be  of  the  order  of  2 x  10^ 

(21  bits)  which  would  require  soma  method,  such  as  gain  switching,  to 
reduce  the  dynamic  range  to  the  manageable  value  of  2'^.  The  system  employs 
uncooled  optics,  but  has  a  cold-shielded  focal  plane  as  well  as  cold  filters. 

Configuration  2  has  the  same  forecptics  and  gross  radiometric  specifi¬ 
cations,  but  is  a  four-port  interferometer  employing  two  dichroics.  The 
diagram  of  Fig.  3.1  shows  that  replacement  of  the  two  plane  mirrors  of  the 
conventional  scanning  Michelson  by  corner  reflectors  results  in  dual  input/ 
output  ports.  Dichroics  in  the  two  output  ports  channel  the  encoded  radia¬ 
tion  to  four  separate  cold-filtered  detector  mosaics,  covering  different 
portions  of  the  total  2.9  -  13  pm  spectral  range.  This  arrangement  allows 
"simultaneous"  detection  of  up  to  four  species,  and  increased  measurement 
time  (greater  detectability)  for  any  single  species. 

In  the  configuration  shown,  the  second  input  port  allows  illumination 
of  the  back  side  of  the  beamsplitter  by  a  reference  blackbody  (although  the 
illuminating  source  could  instead  be  a  different  portion  of  the  external 
scene).  The  effect  is  to  subtract  the  reference  blackbody  spectrum  from 
the  input  spectrum  (or  to  subtract  the  spectra  of  displaced  IFOV's),  which 
corresponds  to  a  potentially  large  reduction  in  the  dynamic  range  of  the 
i nterferogram. 
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CORNER 

REFLECTOR 


In  the  case  where  both  input  ports  receive  external  radiation,  the 
configuration  performs  an  optica!  subtraction  of  the  radiances  reaching  a 
given  detector  from  different  parts  of  the  scene.  Thus,  it  is  not  required 
to  subtract  the  outputs  of  different  detectors  to  suppress  the  background, 
as  in  Configuration  1.  However,  the  noise  in  the  difference  spectrum  (the 
signal  to  be  analyzed)  is  the  same  for  both  configurations:  a  fl  increase 
in  the  NESR  of  Configuration  2  results  from  doubling  the  total  band  photon 
radiance  on  the  BLIP  detector  (the  effect  of  the  second  input  port),  whereas 

in  Conf i gurat ion  1  the  effective  increase  is  the  result  of  subtracting 

k 

two  independent  noise  signals. 

If  the  second  input  port  is  a  blackbody  reference  channel,  and  the 
continuum  portion  of  the  external  scene  spectrum  is  approximately  a  blackbody 
at  the  same  temperature  as  the  internal  reference,  then  it  is,  again,  not 
necessary  to  subtract  the  outputs  of  different  detectors.  The  decision  to 
use  an  internal  versus  external  reference  should  be  based  on  the  expected 
spatial  homogeneity  of  the  background,  the  expected  amount  of  spectral 
structure  in  the  background  (this  determines  the  maximum  reduction  in 
dynamic  range  using  the  blackbody  reference),  the  ability  to  predict  the 
average  spectral  radiance  of  the  background,  and  the  impact  of  the  decision 
on  optical  design  complexity. 

The  MDQ  values  given  in  Section  I  and  Vol .  II  represent  Baseline  Config¬ 
uration  2.  The  assumed  radiometric  properties  (reflectivities,  emi ss i v i t ies , 
etc.)  of  the  optical  elements  of  Configuration  2  are  given  in  Section  1. 

Actually,  the  photon  fiux  from  external  sources  in  Conf i gurat ion  2  would 
be  doubled,  but  the  flux  from  internal  sources  would  increase  by  somewhat 
less  than  a  factor  of  2.  The  assumption  of  v2  increase  in  NESR  is, 
therefore,  somewhat  pessimistic. 


It  should  be  noted  that  the  effects  on  NESR  of  port  4  are  not  included  in 
the  NESR  values  tabulated  in  Section  1  and  Vol .  II;  rather,  the  /2  increase 
is  attributed  to  IFOV  subtraction  in  the  data  processing.  If  the  band 
photon  radiances  (BPR's)  given  in  these  sections  are  used  to  evaluate 
detector  performance  at  different  background  flux  levels  for  Configuration  2, 
the  BPR's  should  be  increased  by  a  factor  of  2.  Note,  also,  that  the  scan 
time  is  fixed  at  the  minimum  value  corresponding  to  background-limited 
performance  (min  t^)  in  order  to  facilitate  scaling  via  the  nomographs  given 
in  Section  1.4. 
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SECTION  4 


EFFECTS  OF  ATMOSPHERIC  TRANSMITTANCE 
(RADIANCE'  VARIATIONS  AND  OFF-NADIR  VIEWING 


4 .  1  Summary 

The  minimum  detectable  amounts  of  the  target  gases  depend,  in  part,  on 
the  transmittance/emission  properties  of  the  intervening  atmosphere,  which, 
in  turn,  are  determined  by  the  atmospheric  composition  and  number  of  air 
masses  in  the  1 i ne-of-s ight .  We  recall  that  the  minimum  detectable  quantity 
(MDQ)  corresponding  to  specified  detection  and  false  detection  probabilities 
is  inversely  proportional  to  a'/  (/2  NESR)  where  o'  is  the  standard  deviation 
of  the  normalized  contrast  spectrum  TvQtgv »  anc*  NESR  is  the  average  system 
noise  equivalent  spectral  radiance  over  the  detection  band.  The  atmospheric 
transmittance  in  the  regions  of  target  emission  lines  influences  o', 
while  the  atmospheric  band  photon  radiance  (and  its  background  absorption) 
affect  the  system  NESR.  In  this  section  we  consider  the  effects  on  MDQ  of 
substituting  the  Midlatitude  Winter  model  atmosphere  for  the  Midlatitude 
Summer  model,  and  the  effect  of  increasing  the  nadir  view  angle  from  0  to 
45  degrees.  Because  of  the  amount  of  computation  involved,  we  evaluated  the 
MDQ  variations  for  only  six  of  the  15  selected  detection  bands.  However, 
these  are  representative  of  the  range  of  variations  to  be  expected  for  all 
detection  bands. 

MDQ  was  defined  as  the  minimum  detectable  target  gas  column  thickness 
per  unit  difference  between  the  average  spectral  radiances  of  blackbodies  at 
the  target  and  background  temperatures.  As  long  as  column  thickness  (number 

J. 

The  choice  of  NO2  detection  band  was  changed  after  this  work  was  performed 
on  the  basis  of  new  information  on  band  strengths.  This  section  treats  the 
Vj  + band,  out  the  preferred  detection  band  is  NO2  (vj)  • 
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of  observed  molecules  per  unit  area)  is  measured  along  the  1 i ne-of-s i ght 
and  the  target  always  fills  the  IFOV,  there  are  no  geometrical  "range- 
squared"  or  "cos6"  effects  in  calculating  the  MDQ.  Of  course,  the  "quantity 
available  for  detection"  can  depend  on  the  viewing  angle  and  range,  depend¬ 
ing  on  the  details  of  the  target  plume  and  viewing  geometries.  Hence, 
minimum  detectable  mass  flow  rates  or  effluent  densities  may  depend  explic¬ 
itly  on  the  nadir  view  angle,  as  well  as  on  the  calculated  MDQ.  For 
example,  if  the  MDQ,  as  defined  above,  happened  to  increase  by  the  factor 
1.1  in  going  zero  to  45  degrees,  and  the  target  gas  were  represented  as  a 
slab  parallel  to  the  earth's  surface,  then,  since  the  quantity  available  for 
detection  increases  by  the  factor  cos  45°  =  /2,  the  minimum  detectable 
volume  concentration  would  change  by  the  factor  1.1//2;  i.e.,  it  would 

actual ly  decrease. 

The  MDQ,  as  defined  above,  varies  as 


MDQ 


NESR 
o  ' 


(J),/2/V 


1/2, 


=  [BPR-0.0228  +  J  1  "Vo' 


(1) 


where  J  is  the  total  photon  flux  density  on  the  detector,  is  the  photon 
flux  density  contribution  of  internal  sources,  and  BPR  is  the  band  photon 
radiance  of  the  external  scene  (atmosphere  plus  background).  The  factor 
0.0228  converts  band  photon  radiance  to  detector  photon  flux  density  for 
the  baseline  system  (op?  cs  transmission  *  0.261  and  F^  =  3-0).  Thus,  the 
MDQ's  calculated  for  the  Midlatitude  Summer  model  atmosphere  and  nadir 

See  Eq .  (23)  in  Section  2.5. 


1 


1 


i 

1 
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viewing  can  be  scaled  easily  to  other  cases  provided  that  BPR  and  o'  are 
determined  for  the  other  cases. 

The  model  atmospheres  most  readily  useable  for  FASCODI  computations  of 
(or  o')  and  in  computations  of  the  band  photon  radiance  are  the  six 
incorporated  into  the  LOWTRAN  series  of  computer  codes.  The  most  signifi¬ 
cant  differences  among  the  models  are  the  total  vertical  water  vapor  amounts, 
which  are  listed  in  Table  4.1  in  units  of  precipitable  centimeters  (pr.  cm). 
Note  that  there  is  a  ten-fold  variation  in  total  1^0  between  the  wettest 
and  dryest  models,  and  that  the  two  models  selected  to  represent  atmospheric 
effects  in  the  MDQ  (the  Midlatitude  Summer  and  Winter  models)  represent 
total  moisture  amounts  of  3.0  and  0.87  pr.  cm.  (which  differ  by  a  factor  of 
3.5).  These  variations  can  be  expected  to  have  significant  effects  on  the 
MDQ.  for  molecules  whose  detection  bands  coincide  with  regions  of 
absorption.  By  comparison,  a  nadir  viewing  angle  change  from  0  to  45  degrees 
causes  a  relative  increase  in  air  mass  of  only  Jl  for  all  atmospheric  species. 

Table  4.1  Vertically  Integrated  Moisture  Content  w 
of  the  LOWTRAN  Model  Atmospheres 


:  1  NO. 

Name 

w  (pr,  1 

1 

T  rop i ca 1 

4.2 

2 

Midlatitude  Summer 

3.0 

4 

Subarctic  Summer 

2.1 

6 

1962  U.S.  Standard 

1  .43 

3 

Midlat i tude  Winter 

0.87 

5 

Subarctic  Winter 

0.42 

The  LOWTRAN  models  encompass  nearly  the  full  range  of  moisture  amounts 
encountered  in  the  atmosphere  and  are  quite  representative  of  different 

q 

latitude  zones.  This  is  illustrated  in  Fig,  4.1,  taken  from  Rosen  et  al  , 
which  shows  the  latitude  variation  of  the  H^O  column  thickness  (after 
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averaging  over  longitude)  for  six  years  of  observation.  Superimposed  on 
Rosen's  plot  are  the  column  thicknesses  represented  by  the  LOWTRAN 
model  atmospheres  (Table  A. I). 


Yt  3  I 


Figure  A. I  Meriodional  Profiles  of  the  Zona!  Mean  Vertically 
Integrated  Moisture  Content  [w]  for  each  of  the 
Years  Studied.  From  Rosen  et  al.® 

The  results  of  the  calculation  of  MDQ  variation  with  model  atmosphere 
and  nadir  viewing  angle  are  summarized  in  Table  4.2.  The  table  gives  the 
lowest  and  highest  values  of  relative  MDQ  which  correspond,  respectively, 
to  the  combinations: 

a. )  Midlatitude  Winter,  nadir  angle  =  zero,  and 

b. )  Midlatitude  Summer,  nadir  angle  =  45°. 

Table  4.2  Maximum  Variation  of  MDQ  for  the  Model 
Atmospheres/Nadir  Angles  Considered 


Molecule , 

Band  (cm 

MDQ  relative 

(lowest) 

to  Ref.  case 

(highest) 

N02 (2850-2935) 

0.75 

1.05 

HCl (2900-3000) 

0.73 

1  .53 

DF(2700-2900) 

0.75 

1 .06 

HBr (2450-2650) 

0.74 

1 .05 

0^(1295-1310) 

0.39 

2.80 

S02( 1090-1210) 

0.53 

1.26 
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The  values  in  the  table  are  relative  to  the  MDQ's  reported  in  Section  1  for 
the  Midlatitude  Summer  case  and  nadir  angle  =  zero;  i.e.,  a  value  of  0.5 
indicates  that  the  MOQ  is  one-half  of  the  previously  reported  value  (the 
molecule  is  "more  easily  detectable"  by  a  factor  of  two). 

It  should  be  noted  that  the  Winter  model  was  used  in  combination  with 
an  assumed  surface  temperature  of  2736K  (f^O  freezing  point):  The  MDQ's 
for  most  of  the  selected  detection  bands  vary  strongly  with  surface  tempera¬ 
ture  (through  the  BPR)  ,  because  most,  of  the  bands  coincide  with  regions  of 
high  atmospheric  transmittance. 

The  results  of  Table  A. 2  show  that  the  MDQ  varies  by  a  factor  of  less 
than  2.4  between  the  lowest  and  highest  values,  except  for  CH^,  where  the 
ratio  of  highest  and  lowest  MDQ  is  approximately  7.  These  relative  varia¬ 
tions  are  explained  by  the  fact  that  the  methane  detection  band  is  in  a 
region  of  strong  water  vapor  absorption,  while  the  other  detection  bands 
contain  only  isolated  1^0  lines  and/or  a  relatively  weak  H20  continuum.  It 
is  recalled  that  CH^  was  one  of  the  molecules  that  is  difficult  to  detect 
through  a  Midlatitude  Summer  atmosphere  (the  minimum  detectable  concentra¬ 
tion  for  a  10-meter  cloud  thickness  and  5°C  temperature  difference  was 
163  ppmV) .  The  implication  of  the  results  in  Table  4,2  is  that  CH^  might 
be  easily  detectable  through  a  dry  Arctic  atmosphere,  and  might  be  impossible 
to  detect  in  useful  quantities  in  the  Tropics.  It  is  possible  that  the 
CH^  Vj  band  in  the  3000-3025  cm  '  region  would  be  a  better  choice  than  the 
band  for  the  tropical  case. 

The  results  shown  in  Table  4.2  for  SO2  represent  an  approximate  upper 
limit  on  the  MOQ.  variation  for  all  molecular  detection  bands  in  the 
8-12.7  um  region  (the  selected  NHj,  HNOj,  anc*  ^3*  bands).  The 


results  obtained  for  NC^,  HC1  ,  DF  and  HBr  are  roughly  indicative  of  the 
results  that  would  be  expected  for  the  ^0,  CO  and  HI  detection  bands. 

These  conclusions  are  based  on  comparisons  of  LOWTRAN  average  transmi ttances 
for  the  different  detection  bands. 

k . 2  Details  of  the  Calculations 

The  six  particular  detection  bands  appearing  in  Table  A. 2  were  selected 
on  the  basis  of  several  considerations.  First,  we  wished  to  represent  mole¬ 
cules  with  sparse  line  spectra  (such  as  HC1  and  DF) ,  as  well  as  those  with 
overlapping  spectral  lines  (such  as  NO2) .  We  also  wanted  to  include  bands 
of  both  high  and  low  atmospheric  transparency,  and  bands  at  short  and  long 
wavelengths  (the  InSb  and  HaCdTe  detector  regions).  To  facilitate  the 
selection  of  bands  we  performed  LOWTRAN  k  calculations  to  obtain,  for  all  of 
the  15  detection  bands,  values  of  average  spectral  transmittance  for  both 
the  Midlatitude  Winter  and  Summer  model  atmospheres  and  for  nadir  view 
angles  of  zero  and  45  degrees.  These  results  are  shown  in  Table  4.3, 

Table  4.3  shows  that  the  average  band  transmi ttances  are  generally 
greater  than  33  percent,  with  the  exceptions  of  the  CO2  band  (which  is  some¬ 
what  of  a  special  case)  and  the  HF  and  CH^  bands,  which  are  in  regions  of 
strong  H2O  absorption.  It  is  interesting  to  note  that  the  average  transmit¬ 
tance  is  higher  over  the  HF  band  than  it  is  over  the  CH^  band  even  though 
we  found  that  CH^  was  detectable  and  HF  was  not.  One  of  the  reasons  for 
this  apparent  inconsistency  is  that  the  average  (low-resolution)  transmission 
is  not  representative  of  the  transmissions  at  the  line  centers  in  the  case  of 
a  molecule  with  widely  separated  lines.  In  fact,  the  spectra'  contrast  plot 
given  in  Vol.  II  for  HF  at  300°K  contains  only  the  PI  1  and  P12  lines  of  the 
1-0  band.  The  detection  band  also  includes  the  P 1 3  and  P)4  lines  at 


3326.21  cm  *  and  3369.90  cm  but  these  apparently  coincide  with  very  strong 
atmospheric  lines,  and,  hence,  do  not  appear  in  the  contrast  spectrum. 

The  LOWTRAN  4  average  transmi ttances ,  although  not  accurate  indicators 
of  available  signal  radiance  for  the  HF,  HC1,  DF,  HBr,  CO  and  HI  bands,  are 
still  useful  for  determining  which  of  these  bands  are  most  likely  to  be 
affected  by  atmospheric  variations.  Based  on  the  values  in  Table  4. 3,  and 
the  other  considerations  mentioned  above,  we  selected  N02,  HC1,  DF,  HBr,  CH^ 
and  S02  as  a  subset  of  detection  bands  that  should  adequately  represent  the 
effects  of  atmospheric  and  nadir  angle  variations. 

The  low-resolution  spectral  transmi ttances  obtained  from  LOWTRAN  4  (in 
the  process  of  computing  the  average  band  transmi ttances)  were  plotted  for 
the  two  extreme  cases:  Midlatitude  Winter,  vertical  viewing  (maximum  spec¬ 
tral  transmittance),  and  Midlatitude  Summer,  45-degree  viewing  (minimum 
transmittance).  The  16  selected  detection  bands  and  the  several  other  candi¬ 
date  bands  considered  are  indicated  in  the  spectra  (Figure  4.2).  This  figure 
is  a  convenient  display  of  the  locations  of  the  detection  bands  relative  to 
the  major  atmospheric  absorption  bands. 

The  relative  MDQ's  for  the  four  different  cases  were  obtained  by  evalu¬ 
ating  Eq.  (1)  and  normalizing  with  respect  to  the  value  obtained  for  the 
reference  case:  Midlatitude  Summer  and  nadir  view  angle  6  =  0.  The  results 
will  be  displayed  in  the  matrix  format  defined  in  Fig.  4.3.  Note  that  two 
values  of  relative  MDQ,  are  given  for  each  of  the  winter  cases.  The  values 
in  parentheses  correspond  to  an  assumed  surface  temperature  of  300°K,  the 
same  value  used  with  the  Midlatitude  Summer  model  atmosphere.  The  other 
values  correspond  to  an  assumed  surface  temperature  of  273°K  (32°F) ,  which 
is  a  more  realistic  surface  temperature  for  a  winter  model.  By  providing 
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gure  4. 


1 

i 


2  Atmospheric  transmittance  of  a  vertical  path  through  the 
Midlatitude  Winter  model  atmosphere  (upper  curve)  and  a  path 
at  45  degrees  zenith  angle  through  the  Midlatitude  Summer 
model  atmosphere  (lower  curve).  Figure  4.2a  covers  the 
600-1200  cm  '  region,  Fig.  4.2b  the  1200-1800  cm  '  region, 
Fig.  4.2c  the  1800-2400  cm  '  region,  Fig.  4. 2d  the  2400- 
3000  cm  '  region,  and  Fig.  4.2e  the  3000-3600  cm  1  region. 
The  selected  detection  bands  (solid  bars)  are  indicated,  as 
well  as  the  other  candidate  detection  bands  considered  for 
N02>  N20,  CH^  and  $02  (dashed  bars).  The  transmittance 
spectra  were  computed  using  the  LOWTRAN  4  code.2* 
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Figure  4.2a  The  600-1200  cm  region 
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Figure  4.2b  The  1200-1800  cin  region 
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Figure  4.3  Format  of  Table  4.4  for  Relative  MDQ  Values 


Midlatitude  Atmosphere 
WINTER  SUMMER 


the  results  for  both  temperatures,  we  can  assess  how  much  of  the  MDQ  varia¬ 
tion  results  from  atmospheric  .ariations  as  compared  to  background  tempera¬ 
ture  variations.  It  is  emphasized  that  background  temperature  variations 
affect  the  MDQ  only  through  the  system  NESR;  the  difference  between  target 
and  background  temperature  is  of  consequence  only  when  MDQ  is  converted  to 
minimum  detectable  column  thickness  or  minimum  detectable  volume 
concentrat i on , 

The  final  results  are  displayed  in  Table  4.4.  In  addition  to  relative 

MDQ,  we  show  the  corresponding  band  photon  radiances  BPR  in  photons/ 

2 

sec  cm  sr  and  the  standard  deviation  o'  of  the  contrast  spectral  radiance; 
these  are  shown  as  pairs  of  matrix  elements  with  the  uppermost  value  equal 
to  tne  BPR.  Beneath  each  BPR,  o'  matrix  we  give  the  value  of  the  photon 
flux  density  on  the  detector  due  to  internal  sources.  These  additional 
values  allow  recalculation  of  the  total  photon  flux  density  on  the  detector, 
the  system  NESR,  and  required  dynamic  range. 

Note  that  for  the  first  four  molecules  in  Table  4.4  the  background 
temperature  has  a  much  larger  effect  on  MDQ  than  the  model  atmosphere.  : 

SC>2  the  model  atmosphere  and  background  temperature  both  have  siun;*  • 
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Table  4.4  Relative  MDQ,  Band  Photon  Radiances,  and  o'  for  Hidiatitude  Winter 
and  Summer  Model  Atmospheres  and  for  0  •  0  and  45  Degrees 


Molecule 
(Detection  Band) 


N02 (2850-2935) 


HC 1(2900-3000) 


DF (2700-2900) 


HBr (2450-2650) 


Relative  MDQ,  Band  Photon  Radiance,  and  o' 


0.75(0.90) 

1.28(3.54)  x  1013 
-20 

1.0  2.22  x  10 

3.41  x  1013 
1.99  x  10"20 

0.77(0.92) 

1,05  1.29(3.41)  x  1013 

3.26  x  1013 

2.17  x  10"20 

1.88  x  10‘20 

8.5  > 

;  ,o- 

,  1.10(2.71)  x  1013 

2.56  x  1013 

0.73(0.84) 

crt 

■ 

o 

X 

00 

o 

• 

9.81  x  10'20 

1.08(1.23) 

1,53  1.11(2.53)  X  JO13 

2.37  x  1013 

7.94  x  lo“20 

6.32  x  10'20 

8.0  x 

lo" 

0.43(1.30)  x  1011* 

1.23  x  lO14 

0.75(0.92) 

1.0  1.04  x  10-19 

9.39  x  10'20 

0.78(0.94) 

1,06  0.42(1.26)  x  10U 

1.18  x  1011* 

1.00  x  10‘19 

8.78  x  10-20 

3.0  > 

1012 

1.05(3.40)  x  lO14 

14 

3.39  x  10IH 

0.74(0.91) 

1.0  2.31  x  »0"20 

2.10  x  10_ZO 

0.76(0.93) 

1,05  1.08(3.32)  x  lO1* 

3.28  x  10‘* 

-20 

-20 

2.25  x  10 

1.99  x  10 

1  O 
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effects,  while  for  CH^,  the  MDQ  Is  determined  primarily  by  the  atmosphere. 
Note  also  that  for  NOj.  DF  and  HBr,  the  variation  of  nadir  view  angle  has 
an  almost  negligible  effect  on  MDQ. 
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SECTION  5 


EFFECTS  OF  INACCURATE  PREDICTIONS  OF 
ATMOSPHERIC  AND/OR  TARGET  GAS  ABSORPTIONS, 

OF  SPECTRAL  INTERFERENCES  AMONG  TARGET  SPECIES, 
AND  OF  INCOMPLETE  BACKGROUND  SUBTRACTION 


In  Section  2,  we  defined  a  procedure  for  obtaining  the  statistics  of 
a  spectral  correlation-detection  process  in  which  the  only  interference  was 
system  (detector)  noise.  That  treatment  ignored  three  "real-world"  effects 
likely  to  degrade  the  system  performance;  these  are 

a. )  imprecise  knowledge  of  the  atmospheric  spectral  transmittance 

and/or  target  gas  spectral  absorption  coefficient  agv, 

b. )  spectral  interferences  (overlap)  between  two  or  more  of  the 

species  to  be  detected,  and 

c. )  different  I ine-of-sight  atmospheric  properties  and/or  back¬ 

grounds  in  the  two  IFOV's  used  to  obtain  the  contrast  spectral 
radiance;  i.e.,  incomplete  subtraction  of  the  net  atmosphere- 
plus-background  spectral  radiance. 

The  analysis  portion  of  Section  2  consisted  of  developing  a  formula 
for  the  coefficient  of  correlation  p  between  the  measured  contrast  spectrum 
and  a  reference  spectrum.  We  found  that  p  could  be  expressed  as  a  ratio 
involving  statistical  descriptors  of  the  reference  spectrum,  the  quant  Ity 
of  target  gas  present,  the  system  NESR,  and  combinations  of  independent 
unit  normal  variates.  Repeated  evaluations  of  this  expression  using 
computer-generated  random  numbers  allowed  us  to  determine  the  statistics  of 
p,  in  particular,  the  detection/false  detection  probabilities  corresponding 
to  various  threshold  levels  on  p,  as  functions  of  the  several  free  param¬ 
eters  in  the  expression.  These  results  were  used  to  determine  minimum 
detectable  target  ges  quantities  corresponding  to  a  baseline  system  n£5W 
and  stated  levels  of  detection  probability  and  false  detection  probability. 
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In  this  section  we  will  redefine  the  contrast  spectral  radiance  and 
reference  spectrum  to  include  the  three  effects  described  above,  and  present 
corresponding  generalized  expressions  for  p.  We  will  then  present  the  results 
of  computations  that  show  the  relative  changes  in  minimum  detectable  quantity 
corresponding  to  assumed  models  of  the  interference  effects.  Derivations  of 
the  generalized  equations  for  p  are  not  included  since  these  follow  the 
same  basic  approach  used  in  Section  2  (and  Appendix  A)  for  the  ideal  case. 

5. 1  Real-World  Effects  in  the  Observed  Contrast  Spectral  Radiance 

Let  N,  with  subscripts  1  or  2,  denote  the  total  upwelling  spectral 
radiance  observed  in  two  different  iFOV's  (subscripts  v  and  g  used  in  pre¬ 
vious  sections  will  be  omitted  to  simplify  the  notation;  e.g.,  tv  will  be 

written  as  t,  tj  or  Tj  and  agv  as  a).  The  first  IFOV  contains  no  target  gas, 

and  observes  spectral  radiance 

N,  -  N*  +  t jN *  ,  (1) 

where  N®  denotes  the  contribution  of  the  atmosphere,  and  TjNj  the  contribu¬ 
tion  of  the  ground  surface  (background)  with  the  effect  of  atmospheric  trans¬ 
mittance  (tj)  included.  The  second  IFOV,  containing  a  cloud  (slab)  of  the 
target  gas  at  ground  level,  observes  spectral  radiance 

N2  "  N2  +  t2(1-t9)B(T9)  +  t9t2N2  ,  (2) 

where  t9  is  the  spectral  transmittance  of  the  gas  slab,  and  B(T9)  is  the 
spectral  radiance  of  a  blackbody  at  the  temperature  T9  of  the  target  gas. 

N®,  N*  and  t2  are  analogous  to  N®,  N*  and  Tj,  and  represent  the  uncontami¬ 
nated  atmosphere  and  background  viewed  by  the  second  IFOV,  which  may  be 
different  than  those  in  the  first  IFOV.  Equations  (1)  and  (2)  correspond  to 
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infinite  resolving  power.  Also,  in  Eq.  (2),  it  is  assumed  that  the  atmos¬ 
pheric  transmittance  between  the  sensor  and  the  top  of  the  cloud  can  be 
approximated  by  the  transmittance  of  the  entire  atmosphere  (12)  when  no  gas 
cloud  is  present;  this  assumption  is  valid  tf  the  physical  cloud  thickness 
is  small  compared  to  one  scale  height  of  the  atmosphere. 

The  contrast  spectral  radiance  AN  is  equal  to  ^"Nj.  If  we  add  and 
subtract  t2N2  to  N2,  perform  the  subtraction  and  collect  terms,  the  result 
is 

AN  -  [N^  +  ^l  “  N,  +  t20-t9)B(T9)  -  ftjN*  -  t9t2N21 
*  N^C'no  gas")  -  Nj  +  TjO't9)  [BCT9)-^] 

2  £  +  t(I-t9)[B(T9)-NS]  ;  «•  =  N2("no  gas")  -N,  .  (3) 

Here,  E  is  defined  as  the  difference  in  net  spectral  radiances  produced  by 
the  (uncontaminated)  atmospheres  and  backgrounds  in  the  two  IFOV's,  and  x 
and  Ns  (note  subscript  "2"  is  omitted)  are  the  transmittance  of  the  atmos¬ 
phere  and  the  intrinsic  background  radiance  in  the  1F0V  containing  the  target. 

If  the  target  gas  slab  is  optically  thin,  Eq.  (3)  reduces  to 

AN  -  £  +  xau[B(T9)  -  NS] 

«  E  +  xau  <B(T9)  -  Ns  >  ,  (*> 

where  <>  denotes  a  spectral  average,  o  is  the  spectral  absorption  coeffi¬ 
cient  of  the  target  gas,  and  u  is  the  optical  thickness  (column  density)  of 
the  target  cloud.  To  the  extent  that  the  spectral  radiance  difference  E  can 
be  tolerated  as  an  "error"  in  the  spectral  measurement,  and  to  (previously 
denoted  Tv<*gv)  can  be  predicted,  the  detectable  quantity  is 
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D  =  u<B(T9)-Ns>  .  (5) 

However,  if  we  acknowledge  that  the  measurement  contains  spectral  noise  n, 
and  that  the  target  may  occupy  a  fraction  F  of  the  IFOV,  the  contrast 
spectra]  radiance  becomes 

AN  *  £  +  toFu<B(T9)  -Ns>  +  n 

*  Oto  +  £  +  n  ,  (6) 

where  the  detectable  quantity  is  now 

D  -  Fu<B(T9)-Ns>  .  (7) 

As  in  the  Ideal  case  (Section  2),  u  can  be  obtained  from  the  inference  of 
0,  provided  F  and  <B(T9)-NS>  can  be  inferred  independently.  If,  in  addi¬ 
tion,  the  physical  thickness  of  the  target  cloud  can  be  estimated,  then  the 
gas  density  or  volume  concentration  can  be  determined. 

Note  that  the  definition  (7)  of  the  detectable  quantity  is  the  same  as 
the  one  given  in  Section  2,  except  that  the  background  is  assumed  to  have 
an  arbitrary  radiance  spectrum  Ns  rather  than  a  blackbody  spectrum  B(TS). 

In  order  to  treat  0  as  a  constant  (spectral  average),  we  must  require  that 
Nj  =  Ns  not  have  spectral  variations  with  structure  comparable  to  that  of  ta 
or  £,  which  is  likely  to  be  the  case.  Of  course,  £  will  have  spectral  struc¬ 
ture  comparable  to  that  of  to,  since  a  component  of  £  is  the  difference  in 
atmospheric  spectral  radiances  in  the  two  IFOV's  (N*-N®) .  It  is  likely  that 
this  component  of  £  will  degrade  detection  performance  much  more  than 
possible  differences  in  background  surface  spectral  radiances. 
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Suppose  there  are  other  (unanticipated)  gas  species  present  in  the 
target  IFOV  that  have  absorption  lines  or  bands  that  overlap  those  of  the 
target  species:  let  Dj ,  D2,  ...  denote  their  detectable  quantities  (Dn  ■ 
Fun<B(T9)  -  Ns  >) ,  and  Op  a, 2 *  •••  their  spectral  absorption  coefficients. 
Then  the  actual  contrast  spectral  radiance  is 

Dt  °2 

AN  *  xD(a  +  a |  -jp  +  02  "jj~  +  . ..)  +  £  +  n  .  (8) 


In  the  spectral  correlation  detection  process  proposed  in  this  study, 
we  attempt  to  detect  the  species  of  interest  by  thresholding  on  the  normal¬ 
ized  correlation  between  the  measured  contrast  AN  and  a  reference  spectrum 
r'a',  where  r'  and  a'  represent  estimates  of  t  and  a.  Specifically,  the 
detection  process  consists  of  thresholding  on  the  value  of 


covariance  of  AN  and  j'a' 


(std  5  standard  deviation) 


It  is  evident  that  the  detection  probability  and  false  target  rejection 
capability  will  both  be  diminished  by  the  noise  term  n,  by  the  difference 
between  xa  and  x'a',  by  the  presence  of  interfering  species  (the  terms 
cxj  Dj/D,  c»2  D2/O,  ...  in  Eq.  (8)),  and  by  the  residual  atmosphere-background 
spectral  radiance  £.  Equations  8  and  9  may  be  used  to  re-evaluate  the 
statistics  of  p  when  one  or  more  of  the  interference  effects  are  included, 
and  then  to  determine  the  corresponding  increase  in  minimum  detectable 
quantity  for  fixed  detection/false  detection  probabilities.  The  following 
subsection  concerns  the  development  of  parameterized  expressions  for  p; 
l.e.,  forms  which  do  not  require  computer  simulations  involving  the  actual 
spectra  xa  and  x'a". 
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5.2  Generalized  Expressions  for  the  Coefficient  of  Correlation  p 

Consider  that  the  interference  effects  referred  to  above  as  a,  b,  and  c 
each  have  "true"  and  "false"  states,  which  we  can  denote  as  I  and  0, 
respectively.  The  number  of  different  combinations  of  states  for  ail  three 
effects  is  therefore  ■  8.  We  can  label  the  combinations  by  the  numbers 
0  through  7.  and  readily  identify  the  specific  states  by  the  binary  equi¬ 
valents  of  the  labels: 


Combination 

0 

1 

2 

3 

4 

5 

6 
7 


Effects 

a  b  c  (1  -  "true",  0  -  "false") 

000  a (true)  means 
001  to  estimated  by  x'a' 

0  1  0 

0  1  1  b(true)  means 

D1 

1  0  0  a  -►  a  +  a,  -jj-  +  ... 

1  0  1 

1  1  0  c(true)  means 

111  C  *  0 


We  derived  new  expressions  for  p  for  combination  7  (all  effects  present 
simultaneously)  and  for  combinations  1,  2  and  4  (each  effect  acting 
independently).  The  expressions  for  the  other  combinations  were  obtained  by 
inference.  As  stated  earlier,  we  will  present  only  the  results  of  the 
derivations. 

It  is  convenient  to  first  present  the  results  obtained  for  combination 
7,  since  all  other  results  can  be  obtained  by  simple  elimination  of  param¬ 
eters  and/or  by  substitution  of  degenerate  forms  for  the  equations  defining 
parameters  in  the  expressions  for  p.  It  should  be  remembered  that  p 
generally  assumes  different  forms  for  D  ■  0  (the  one  used  to  compute  false 
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detection  probability)  and  for  D  i1  0  (the  one  used  to  compute  the  detection 
probabf 1 ity) . 

5.2.1  Results  for  Combination  7 

For  this  combination  we  find  that  the  coefficient  of  spectral  correla¬ 
tion  between  AN  and  t 'ct'  is 


(10) 


where  the  x.  are  independent  unit  normal  random  numbers,  M  is  the  number  of 
resolved  spectral  elements  comprising  AN  and  T'a',  and  where 


P 


z 


!°i  *  °u 


■0{  * 


> 


ri[°o2+  2aC 
>[o02+  2o0c  +  of^/(S2  NESS) 


(11) 


(12) 


2  2 

oQ  »  <(tB-u)  >  *  variance  of  t0  (actual  ta-type  spectrum)  (13) 

2  2 

<Tj  =  <(r'o'-u')  >  *  variance  of  r'a'  (reference  to  spectrum)  (14) 

2  2 

a-  ■  <(£/D-p  )  >  -  variance  of  C/D  (normalized  spectrum  error  in  (15) 
*  *  background  subtraction) 

°01  “  <(tB-u)  (x'a'-u')>  *  covariance  of  t6  and  t'a'  (16) 

°0£  "  <(x6"u) (5/0-u^)>  "  covariance  of  x6  and  £/D  (17) 


*  <(t'o'-jj')  (5/0-y^)>  •  covariance  of  r'a'  and  £/D. 


08) 


For  0  +  0,  6  is  defined  by 


6  *  a  +  (Dj/D)aj  +  (D2/D)c»2  +  •••  (19) 

while  for  "0  »  0" 

0  »  (Dj/D)aj  +  (Dj/Djc^  +  ...  ;  D  =  any  value  (19a) 

In  the  case  D  *  0,  which  is  used  to  obtain  the  false  detection  proba¬ 
bility,  p  (and  also  j5  and  z)  are  independent  of  D  (see  Eqs.  10-12,  15,  17, 
18,  19a).  Thus,  any  value  of  D  may  be  used  in  Eqs.  12,  15,  17,  18  and  19a. 

In  the  definitions  (13)“(l8),  the  symbols  y,  y'  and  y^  denote  the 
means  of  the  spectral  quantities  included  in  the  same  set  of  parentheses. 

In  particular,  y  is  the  mean  of  tg,  and,  hence,  it  assumes  different  values 
for  the  cases  D  )<  0  and  "0  -  0". 

The  parameter  j5 ,  like  p,  has  the  form  of  a  coefficient  of  correlation. 
In  fact,  we  see  from  Eq.  (10)  that  it  is  the  value  assumed  by  p  when  there 
is  no  system  noise  (all  x.  ■  0) .  Likewise,  we  will  find  that  $  ■  1  when 
effects  a,  b  and  c  are  absent. 

It  is  convenient  to  consider  next  the  combinations  6,  3  and  2,  which, 
like  7,  have  target  species  overlap  as  one  of  the  interference  effects. 

5-2.2  Results  for  Combination  6 

This  combination  is  the  same  as  7  except  that  c  ■  0.  That  is,  £  *  0, 
and  hence  any  variances  or  covariances  with  subscript  £  are  equal  to  zero. 
Specifically,  Eqs.  (II)  and  (12)  reduce  to 


°0!/(o1°0) 

(11a) 

0o0/(^2  nBE) 

(12a) 

100 


2  2 

Oq  ,  Oj  and  cQ1  are  defined  by  Eqs.  (13),  (14)  and  (J6).  Equations  (19) 
and  (20)  are  still  applicable. 

5.2.3  Results  for  Combination  3 

This  combination  is  the  same  as  7  except  a  =  0.  That  is,  t'a'  +  tq  in 
Eqs.  (14),  (16)  and  (18).  The  other  equations  are  unchanged. 

5.2.4  Results  for  Combination  2 

Combination  2  is  the  same  as  6  except  a  =  0.  That  is,  Eqs.  (11)  and  (12) 
reduce  to  (11a)  and  (12a),  and  r'a'  -*■  to  in  Eqs.  (14),  (16)  and  (18). 

We  next  consider  combinations  5  and  1,  which  are  distinguished  by  the 
property  c  ■  1  and  b  =  0  (different  atmospheric  and/or  background  radiances 
in  the  two  IFOV's,  but  no  spectral  interferences  among  target  species). 

Note  that  when  a  particular  detection  band  contains  no  overlapping  target 
species  (or  when  overlapping  species  are  not  expected  to  occur 
simultaneously),  the  only  combinations  of  interest  are  5,  1,  4  and  0. 

5.2.5  Results  for  Combination  5 

Since  there  are  no  spectral  interferences,  ctj  =  a2  =  ...  =0,  and  there¬ 
fore  Eq.  (19)  gives  6  ■  a  for  D  ^  0  and  6=0  for  D  =  0.  Alternately,  we  may 
write  the  entire  set  of  equations  in  the  form 


(20)  (same  as  10) 
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where,  if  0  +  0 


5-2.6  Results  for  Combination  4 
When  D  =  0,  we  find  that 


P 


xi 


(31) 


For  D  t  0,  Eq.  (20)  or  (10)  is  applicable.  Combination  4  is  the  same  as  5 
except  c  -  0  (£  =  0) ;  hence,  from  Eqs.  (21)  and  (22),  we  obtain 


(5  =  o0]/(a0a,) 

(32) 

z  -  aQD/(/2  NESR) 

(33) 

as  the  parameters  to  use  in  Eq.  (20).  oQ2,  a2  and  oQ1  are  given  by 
Eqs.  (25),  (26)  and  (28). 


5*2.7  Results  for  Combination  1 

This  combination  is  the  same  as  5  except  that  a  *  0.  That  is,  r'a'  -*■ 
xa  in  Eqs.  (26),  (28)  and  (30).  With  this  substitution  aQ2,  a,2  and 
become  identical,  as  do  aQ^  and  o^.  Thus,  Eqs.  (21)  and  (22)  can  be 
written 


70  c 


°Q  +  2a0S  + 


z 


°o2  + 


2o 


0? 


+  o. 


(34) 


(35)  (same  as  22) , 


which  apply  when  0  ^  0.  Equations  (23)  and  (24)  are  applicable  when  D  =  0, 


103 


except  that  an^  0|  °q: 


V(v5> 


(36) 


kz 


TT 


z  =  o^D/ (/?  NESR) 


(37) 


5.2.8  Results  for  Combination  0 

This  combination  is  the  same  as  4  except  that  a  *  0.  That  is,  i'a'  -*•  ra 
in  Eqs .  (26),  (28)  and  (30).  Then  of  -  oQ)  -  a,2,  and  Eqs.  (32)  and  (33) 
reduce  to 


p  =  1 


z  -  o00/(vTNESR)  . 


(38) 

(39) 


Equations  (38),  (39)  and  (20)  define  p  for  D  ^  0.  For  D  *  0,  p  is  given  by 
Eq.  (31).  These  are  exactly  the  results  obtained  in  Section  2  for  the 
case  a  *  b  *  c  =  0  (note  that  oQ  *  Oj  was  previously  denoted  by  a'). 

5.3  General  Influence  of  Real -World  Effects  in  p.  in  the  Detection/False 

Detection  Probabilities,  and  in  the  Minimum  Detectable  Quantities 

To  understand  how  the  "real  world"  effects  influence  the  MDQ's,  we 
must  first  consider  how  the  probability  density  function  (pdf)  of  the 
absolute  correlation  coefficient  |p|  behaves  with  detector  noise  and  with 
the  types  of  interference  described  above.  We  are  able  to  describe  this 
behavior  in  general  terms  with  the  aid  of  hindsight  gained  from  realistic 
experiments  on  the  computer,  and  to  show  that  one  of  the  real-world  effects 
will  essentially  eliminate  any  possibility  of  detecting  some  of  the  species 
of  interest,  unless  the  detection  scheme  described  in  Section  2  is  modified. 
However,  the  required  modification  is  very  simple  and  appears  to  be  quite 
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effective;  that  is,  it  results  in  only  a  slightly  increased  MDQ  for  the  sev¬ 
eral  species  considered.  Section  5.4  describes  the  modified  detection  proce¬ 
dure,  and  Section  5.6  presents  some  revised  MDQ's  for  a  worst-case  situation. 
Section  5.5  describes  the  interference  model  used  to  compute  the  revised  MDQ's, 


Figure  5.1  Pdf's  of  jp|  for  the  Ideal  Case 


Figure  5.1  shows  how  the  pdf  of  |p|  behaves  with  signal-to-noise  z  in 
the  case  when  detector  noise  is  the  only  interference  (the  "ideal  case",  or 
case  0,  which  was  treated  in  Section  2;  the  detection  signal-to-noise  z 
was  denoted  Z  in  that  section).  The  dashed  curve  in  the  first  part  of  the 
figure  is  the  pdf  when  no  gas  is  present  (D*0),  and  the  solid  curve  is  the 
pdf  when  Dj<0.  Note  that  this  is  a  schematic  representation;  the  solid 
curve  will  generally  not  be  symmetrical  about  its  modal  value.  The  area 
under  each  curve  is  unity.  The  area  under  the  solid  curve  at  |p|  greater 
than  some  threshold  Rj  is  the  detection  probability,  while  the  area  under 
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the  dashed  curve  at  |p|  > Rj  is  the  false  detection  probability.  The  first 
part  of  the  figure  corresponds  to  low  signal-to-noise  z  or  to  0  ■  the  MDQ. 
For  example,  if  M*1000  resolution  elements,  a  z-value  of  0.13  gives 
detection/false  detection  probabilities  of  approximately  0.95  and  0.01, 
respectively,  when  the  threshold  Rj  is  0.075.  Generally,  as  z  is  increased 
the  solid  curve  moves  towards  higher  |p|,  and  becomes  narrower  (of  course, 
it  maintains  unit  area).  The  dashed  curve,  however,  is  unchanged  since  it 
is  independent  of  z  (see  Eq.  (31)1.  For  high  z-values  we  can  choose  a 
threshold  that  gives  essentially  unity  detection  probability  and  zero 
false  detection  probability.  For  infinite  z  the  solid  curve  becomes  a  unit 
delta  function  at  jp|  =1.  If  the  detector  noise  is  zero  the  dashed  curve 
becomes  a  delta  function  at  |p|  =0. 

Figure  5.2  compares  the  pdf's  for  the  ideal  and  non-ideal  cases,  where 
"non- ideal"  means  that  one  or  more  of  the  real-world  effects  can  be  present. 
In  general,  in  the  non-ideal  case,  the  pdf  for  D  =  0  (the  dashed  curve)  and 
the  pdf  for  01*0  (the  solid  curve)  are  associated  with  different  z-values 
and  with  different  values  of  the  parameter  p.  (We  recall  that  p  is  the 
value  of  the  correlation  coefficient  when  there  is  no  detector  noise;  it 
is  unity  in  the  ideal  case  when  0  0) .  The  position  of  the  maximum  in  the 

dashed  curve  depends  on  both  z  and  p,  as  does  the  "width"  of  the  curve; 
the  same  is  true  of  the  solid  curve.  Part  (b)  of  Fig.  5.2  is  an  example  of 
the  appearance  of  the  pdf  of  |p|  for  the  non-ideal  case.  Note  that  both  the 
dashed  and  solid  curves  (in  this  example)  are  shifted  towards  higher  jp|  and 
are  narrower  realtive  to  those  for  the  ideal  case.  The  shift  occurs  when 
the  interferences  in  the  spectrum  have  "accidental"  correlation  with  the 
reference  (t'o')  spectrum.  The  narrowing  of  the  curves  occurs  when  the 
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interference  increeses  the  radiance  observed  in  spectral  resolution  elements 
where  the  gas  of  interest  emits  (or  increases  the  absorption  when  the  gas  of 
interest  is  observed  in  absorption);  i.e.,  when  the  effective  "signal-to- 
noise"  is  increased.  However,  the  situation  depicted  in  part  (c)  of  Fig.  5.2 
can  also  occur.  Here  the  maximum  of  the  dashed  curve  corresponds  to  larger 
( p |  than  the  maximum  of  the  solid  curve,  which  means  that  the  presence  of 
the  gas  of  intc.ast  tends  to  decrease  the  value  of  the  correlation 
coefficient . 

Assume  for  the  moment  that  part  (b)  of  Fig.  5.2  is  the  only  type  of 
non-ideal  situation  than  can  exist,  in  this  situation  it  is  sometimes  possi¬ 
ble  to  choose  a  new  threshold  R2  that  gives  detection/false  detection  proba¬ 
bilities  on  the  order  of  0.95/0.01  even  for  D  equal  to  the  MDQ  of  the  ideal 
case.  However,  in  order  to  achieve  a  0.95  detection  probability  in  the  ideal 
case  using  this  new  threshold,  D  would  obviously  have  to  be  larger  than  the 
original  MDQ;  that  is,  the  pdf  of  |pj  would  have  to  look  like  part  (d)  of 
Fig.  5.2.  For  this  value  of  0  the  pdf  of  |pj  for  the  non-ideal  case  would 
look  like  part  (e)  rather  than  part  (b)  of  the  figure.  The  new  value  of  D 
would  be  the  new  MDQ  and  R2  would  be  the  new  required  threshold  on  |p|. 

The  pdf's  shown  In  parts  (b)  and  (e)  of  Fig.  5.2  (and  also  part  (c)) 
are  the  type  that  one  encounters  in  modelling  a  particular  or  perhaps  worst- 
case  non-ideal  situation.  In  practice,  the  interference  can  range  in  magni¬ 
tude  from  zero  (the  ideal  case)  to  the  worst  case  situation,  according  to 
some  other  probability  density  function.  The  effect  in  part  (b)  is  to 
broaden  the  two  curves  and  shift  the  maxima  towards  smaller  |p|.  The  same 
effect  occurs  in  part  (e) ,  which  tends  to  place  more  of  the  pdf  areas  in  the 
vacant  regions  (0,|p|a)  and  (R2,|p|b);  however,  the  threshold  1 1  still 

result  in  adequate  detection  probability  and  adequate  false  target  rejection. 


Unfortunately,  the  two  non-ideal  situations  depicted  in  parts  (b)  and 
(c)  of  Fig,  5.2  can  occur  with  roughly  equal  (50  percent)  probabilities  for 
the  type  of  interference  represented  by  different  atmospheric/background 
properties  in  the  two  detection  IFOV's.  This  means  that  the  best  performance 
that  can  be  achieved  is  roughly  50  percent  detection  probability  and  50  per¬ 
cent  false  detection  probability,  which,  of  course,  is  an  unacceptable  level 
of  performance. 

5.1*  Modified  Detection  Procedure 

We  will  first  illustrate  how  the  different  situations  displayed  in 
parts  (b)  and  (c)  of  Fig.  5.2  can  occur  when  the  two  IFOV's  contain  different 
amounts  of  water  vapor,  and  then  show  how  a  simple  modification  of  the  corre¬ 
lation  detection  process  can  all  but  eliminate  the  effects  of  this  type  of 
interference. 

Suppose  that  IF0V1  contains  more  1^0  than  IF0V2  (we  recall  that  IF0V2 
contains  the  target  gas).  All  of  the  selected  detection  bands  are  located  in 
either  the  3“5  pm  region  or  the  8- 1 3  pm  region,  where  atmospheric  1^0  has 
weak  continuum  absorption.  Under  normal  conditions  the  tropospheric  tempera¬ 
ture  lapse  rate  is  positive  (temperature  decreases  with  altitude)  and,  hence, 
the  H20  acts  as  a  net  absorber  rather  than  as  an  emitter.  Consequently,  the 
observed  continuum  radiance  level  will  decrease  as  the  1^0  column  thickness 
in  the  IF0V  is  increased.  The  observed  spectral  radiances  for  the  two  IFOV's 
will  look  more  or  less  as  indicated  in  Fig.  5.3a,  assuming  the  target  gas  is 
observed  in  emission  (  is  hotter  than  the  surface). 

The  detection  bands  will  usually  contain  some  strong  atmospheric  spec¬ 
tral  lines,  usually  due  to  1^0.  In  lines  that  are  extremely  strong  ("opaque") 
the  observed  spectral  radiance  will  be  very  small  for  both  IFOV's,  as 
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^  PSm — ' — Target  Gas  Emission 


Spectrum  (b)  adjusted 
to  zero  mean 


Spectrum  (c)  adjusted 
to  zero  mean 


Reference  spectrum 
adjusted  to  zero  mean 

figure  5.3  Effects  of  Strong  Atmospheric  Lines  on  the 
Mean  of  the  Contrast  Spectral  Radiance 
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Indicated  In  Fig.  5.3a.  The  contrast  spectrai  radiance  N2~Nj  will  resemble 
part  (b)  of  Fig.  5.3.  Clearly,  if  the  excess  HjO  is  transferred  from  IF0V1 
to  IF0V2,  the  contrast  spectrum  would  have  the  appearance  of  part  (c)  of  the 
figure.  What  we  wish  to  demonstrate  is  that  parts  (b)  and  (c)  are  the 
respective  counterparts  of  the  pdf's  represented  in  parts  (b)  and  (c)  of 
Fig.  5-2. 

We  note  that  contrast  spectrum  5.3b  has  a  positive  mean,  and  spectrum  5-3c 
has  a  negative  mean.  The  two  spectra  adjusted  to  zero  mean  (an  operation 
prescribed  by  the  spec t  'll  correlation  process)  would  look  like  Figs.  5.3d 
and  5.3e,  respectively.  The  zero-mean  spectra,  of  course,  would  be  identical 
if  the  contributions  o'  strong  atmospheric  lines  were  absent.  In  the  corre¬ 
lation  processing  -rtctrum  (d)  or  (e)  is  correlated  with  the  zero-mean 
reference  spectrum  pictured  in  Fig.  5.3f.  Generally  (d)  and  (f)  will  have 
imperfect  Lit  high  correlation  while  (e)  and  (f)  will  have  Imperfect  but 
relatively  low  correlation.  Low  correlation  is  obtained  in  the  latter  case 
because  small  and/or  negative  values  of  apparent  target  spectral  radiance 
coincide  with  high  spectral  radiance  values  in  the  reference  spectrum 
(Multiplying  the  zero-mean  contrast  radiance  by  -1  would  have  no  effect 
since  we  test  on  the  absolute  value  of  p;  it  is  necessary  to  test  on  |p| 
because  the  target  can  be  observed  in  either  emission  or  absorption).  We 
see  that  the  correlation  would  actually  be  greater  if  IF0V2  contained  no 
target  gas.  This  results  in  the  pdfs  shown  in  part  (c)  of  Fig.  5.2.  When 
the  IFOV  containing  the  emitting  target  has  less  HjO  than  the  other  IFOV  the 
resulting  pdf's  look  like  Fig.  5.2a.  Note  that  correlation  coefficients 
higher  than  those  of  the  ideal  case  can  be  obtained  merely  because  the 
interference  results  in  higher  correlation  signal-to-noise  z. 


An  obvious  solution  to  this  dilemma  is  to  exclude  from  the  correlation 
process  any  spectral  elements  containing  very  strong  atmospheric  absorption 
lines.  That  is,  we  should  assign  zero  values  to  their  observed  radiances 
and  reduce  the  total  number  of  spectral  elements  M  in  proportion  to  the 
number  of  assigned  zeros.  Note  that  the  reference  spectrum  t'o'  is  effec¬ 
tively  zero  in  the  regions  of  opaque  lines;  hence,  the  real  effect  of  modi¬ 
fying  the  correlation  procedure  is  to  prevent  strong  atmospheric  lines  from 
affecting  the  means  of  the  observed  radiance  spectra  for  the  two  IFOVs. 
Clearly,  the  correlation  processing  must  include  subtraction  of  the  mean 
spectral  radiance:  otherwise,  the  correlation  would  be  reduced  by  different 
1^0  continuum  levels  and/or  different  surface  temperatures  in  the  two  IFOVs. 

5-5  Interference  Model  Used  to  Estimate  Revised  MDQ's 

The  sources  of  the  target  gas  emissions  are  not  likely  to  emit,  simul¬ 
taneously,  species  that  have  overlapping  detection  bands,  at  least  for  the 
15  species  considered.  That  is,  interference  type  "b"  need  not  be  considered. 
In  modelling  the  effects  of  interference  type  "c",  we  assumed  that  the  only 
differences  between  the  atmospheres/backgrounds  observed  by  the  two  IFOVs 
are  the  HjO  column  thicknesses  and/or  the  background  temperatures.  In  model¬ 
ling  interference  type  "a",  we  assumed  that  the  estimated  reference  spectrum 
x'a'  corresponds  to  an  incorrect  estimate  of  the  1^0  column  thickness;  i.e., 
that  cT»a,  but  x'tx  only  because  of  the  wrong  1^0  column  thickness. 

One  would  not  expect  interference  type  "a"  to  degrade  system  perform¬ 
ance  nearly  as  much  as  type  "b",  since  the  latter  can  introduce  lines  in  the 
contrast  radiance  spectrum  whose  positions  are  random  with  respect  to  the 
target  gas  lines.  Also,  the  effect  of  different  surface  temperatures  in  the 
two  IFOVs  can  be  greatly  reduced  for  many  of  the  detection  bands  by 
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"pre-whitening"  the  observed  contrast  radiance  spectrum;  that  is,  by 
removing  any  slope  or  curvature  in  its  envelope.  Thus,  in  the  modelling  of 
interference  effects  this  study  emphasizes  the  case  of  different  H^O  amounts 
in  the  two  IFOV's,  Spectral  variations  in  earth  background  emissivity  were 
not  considered,  because  they  probably  have  only  low  resolution  features  that 
can  also  be  eliminated  by  pre-whitening. 

We  do  not  know  of  any  data  on  HjO  vertical  column  thickness  variations 
over  horizontal  distances  of  the  order  of  10  meters.  However,  a  reasonable 
estimate  of  the  maximum  variation  can  be  obtained  by  considering  the  mech¬ 
anisms  involved  in  the  formation  of  fair  weather  cumulus  clouds  (this  point 
of  view  was  suggested  to  us  by  Dr.  J.  Bunting  of  AFGL) .  Generally,  these  low 
clouds  are  formed  when  ground  level  moisture  rises  and  cools  to  the  point 
where  saturation  occurs.  A  reasonable  model  for  the  vertical  H20  density 
profile  (or  dew  point  temperature  profile)  under  the  cloud  is  therefore 
obtained  by  assuming  constant  dew-point  temperature  from  ground  level  up  to 
the  altitude  where  the  actual  temperature  has  lapsed  to  this  value.  For  the 
LOWTRAN  Midlatitude  Summer,  Subarctic  Summer  and  Tropical  model  atmospheres, 
we  find  that  this  altitude  is  very  close  to  one  kilometer,  which  is  the 
height  usually  observed  for  low  cumulus  cloud  bases  in  these  types  of  climate 
(they  can  be  somewhat  higher  for  extremely  dry  climates). 

The  assumption  is  made  that  the  maximum  difference  in  total  H^O  in  the 
two  detection  IFOV's  is  equal  to  the  excess  HjO  below  the  cumulus  cloud,  as 
defined  by  the  above  argument.  That  i s ,  we  assume  that  the  worst  case  con¬ 
dition  is  when  a  cloud  is  about  to  be  formed  at  one  km  in  one  of  the  IFOV's. 
The  calculation  of  the  H20  difference,  using  LOWTRAN  model  atmospheres,  is 
then  quite  simple:  we  merely  determine  the  H20  amount  in  the  0-1  km  layer 


assuming  the  LOWTRAN  1^0  density  varies  linearly  between  0  and  l  km,  and 
subtract  thir.  from  I  he  amount  obtained  nssuminq  the  density  is  equal  to  the 
ground  level  value  between  0  and  1  km.  By  this  procedure  we  find  that  the 
differences  are 

0.15  pr  cm  for  the  Subarctic  Summer  model 

(which  has  a  total  of  2.1  pr  cm) 

0.20  pr  cm  for  the  Midlatitude  Summer  model 
(which  has  a  total  of  3.0  pr  cm) 

0.30  pr  cm  for  the  Tropical  model 

(which  has  a  total  of  4.2  pr  cm) 

In  each  case  the  difference  is  approximately  seven  percent  of  the  total 
column  thickness.  In  modelling  the  effects  o*  the  difference  we  used  the 
Midlatitude  Summer  model  and  a  difference  of  0.2  pr  cm.  One  calculation  was 
made  of  the  effects  of  simultaneous  interferences  "a"  and  "c";  in  this  case 
the  "estimated"  spectral  transmi ttance  represented  an  error  in  estimated  1^0 
amount  equal  to  0.4  pr  cm. 

5.6  Some  Revised  MDQ's 

We  examined  the  effects  of  the  interferences  for  the  preferred  SO2,  HCl 
and  HBr  detection  bands.  We  picked  the  SOj  and  HCl  bands  because  their  MDQ.'s 
are  more  strongly  affected  by  HjO  variations  and/or  nadir  view  angle  than 
the  other  bands,  provided  we  exclude  CH^,  which  is  only  marginally  detectable 
(see  Table  4.2,  Section  4).  Also,  SOj  and  HCl  are  different  In  that  the 
former  has  significant  overlap  between  its  lines,  while  the  latter  has  very 
widely  separated  lines.  We  selected  HBr  as  the  third  band  because  it  con¬ 
tains  very  little  H^O  absorption. 

Interference  effects  in  SOj  detection  were  partially  evaluated  for  the 
original  method  of  data  processing,  but  a  revised  MDQ.  value  was  computed  only 
for  the  revised  method,  which  masks  out  strong  atmospheric  lines.  For  the 


other  two  molecules  we  evaluated  only  the  revised  method.  The  revised 
method  was  implemented  using  the  following  simplified  procedure:  we  first 
examined  the  previously  mentioned  spectral  transmittance  plots  for  the 
Midlatitude  Summer  model  (see  Section  1  and  Vol  .  II),  and  determined  the 
average  value  of  the  spectral  transmittance  envelope.  We  then  selected  a 
somewhat  lower  transmittance  value  to  use  as  a  threshold  in  the  correlation 
processing.  That  is,  spectral  elements  having  predicted  transmi ttances  less 
than  the  threshold  value  were  excluded  when  we  evaluated  the  spectrum 
variances  and  covariances  that  determine  the  correlation  coefficient  p. 

For  example,  the  Vj  band  of  SO2  has  a  spectral  transmittance  envelope  equal 
to  approximately  O.58;  the  selected  transmittance  threshold  is  0,4,  An 
examination  of  the  spectral  transmittance  plot  for  this  band  shows  this 
threshold  will  eliminate  the  effects  of  a  large  number  of  the  atmospheric 
absorption  lines,  except  perhaps  in  the  far  wings  of  the  strongest  lines. 

5.6.1  The  V|  Band  of  SO^ 

We  first  computed  the  normalized  contrast  spectral  radiance 

xa + £/D  without  spectral  masking  to  demonstrate  that  strong  i^O  lines  have 

large  contributions  to  the  contrast.  Figure  5.4a  shows  the  normalized 

spectral  contrast  when  IF0V2  contains  SO2  in  the  quantity  2.66  x  101'  (molec/ 

2  2  "I 

cm  ) (W/cm  srcm  ),  which  is  equal  to  the  MDQ  for  a  dwell  time  of  10  seconds, 
and  IF0V1  contains  0.2  pr  cm  more  than  IF0V2.  Figure  5.4b  is  the  norma¬ 
lized  reference  contrast  spectrum  xa  =  x'a'  plotted  to  the  same  scale  as 
Fig.  5.4a.  In  the  original  method  of  data  processing,  detection  of  SO2  would 
be  based  on  the  correlation  of  parts  (a)  and  (b)  of  the  figure.  Note  that  the 
normalized  spectral  contrast  xa+^/0  will  tend  to  be  independent  of  D  at 
frequencies  where  SO2  is  the  only  emitter,  and  varies  as  1/D  where  f^O 
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Figure  5.*t  Part  (a):  Normalized  contrast  spectral  radiance  tct  +  C/D  for 

1 1  2 

the  V|  band  of  SO^  for  0  =  2.66  x  10  (molec/cm  ) 
(W/cm2  sr  cm  ').  IF0V1  contains  0.2  pr  cm  more  H20 
than  IF0V2. 

Part  (b) :  The  normalized  spectral  contrast  to  plotted  to  the 
same  scale  as  part  (a). 


6 


or  some  other  non-target  species  is  the  only  emitter.  However,  on  the 
normalized  "zero-to-one"  scales  of  the  figures,  the  S02  features  would  tend 
to  become  more  prominent  as  D  increases.  Note,  also,  that  the  H20  lines 
in  the  S02  detection  band  have  a  significant  effect  on  the  mean  value  of 
the  normalized  spectral  radiance  (when  0  equals  the  MDQ.  for  t^  =10  seconds). 

Figures  5.5a  and  5-5b  show  the  corresponding  results  when  the  excess 
0.2  pr  cm  H20  is  in  IF0V2  rather  than  IFOVl,  Now,  the  actual  spectral  con¬ 
trast  is  negative  at  all  frequencies  (the  horizontal  lines  denote  zero 
contrast;  again,  the  scales  are  identical).  It  may  be  noted  that  the 
reference  spectrum  is  somewhat  smaller  than  for  the  case  shown  in  Fig.  5. A. 
This  is  because  the  predicted  reference  spectrum  is  assumed  to  have  no  error 
(x"a^=To),  and  because  t  is  now  smaller  as  a  result  of  the  increased  H20 
above  the  S02  in  IF0V2. 

The  computed  parameters  that  determine  the  statistics  of  |p|  for  the 
two  cases  are: 


FIGURE  5.4 
(Aw  =  0.2  in  1F0V1) 

FIGURE  5.5 
(&w  =  0.2  in  IF0V'} 

M 

1200 

1200 

p(D  *0) 

0.573 

0.110 

z-NESR(D  t  0) 

5.02  x  lo'9 

4. 14  x  10"9 

p(D  =  0) 

0.273 

-0.273 

z-NESR(0  =  0) 

4.28  x  lo'9 

4.28  x  10"9 

We  note  that  the  presence  of  S02  tends  to  increase  |p|  (and,  hence  |p|)  as 
well  as  the  correlation  s ignal -to-noi se  z  in  the  first  case,  while  it  has 
the  opposite  effect  in  the  second  case.  If  we  were  to  use  the  above 
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parameters  to  compute  the  pdf's  of  jp),  we  would  find  that  they  resemble 
Figs.  5.2b  and  5.2c,  respectively. 


Figures  5.6  and  5.7  show  the  result,  in  the  contrast  spectrum,  of 
"masking  out"  spectral  elements  with  transmi ttances  less  than  the  threshold 
transmittance  value,  0.4.  In  these  figures  the  excluded  spectral  elements 
have  been  assigned  a  contrast  value  equal  to  the  mean  contrast  of  the  remain¬ 
ing  elements  (in  the  correlation  processing  they  are  simply  excluded). 

Figure  5.6  represents  the  case  of  0.2  pr  cm  excess  1^0  in  1 FOV I ;  the  number 
of  elements  remaining  after  applying  the  exclusion  criterion  is  926  compared 
to  the  original  1200.  Figure  5.7  represents  the  case  of  0.2  pr  cm  excess 
H2O  in  IF0V2;  the  new  M  is  870.  In  the  revised  correlation  procedure,  the 
926  surviving  elements  of  Fig.  5.6  are  correlated  with  corresponding  elements 
of  Fig.  5.4b;  similarly,  the  870  elements  of  Fig.  5.6  are  correlated  with 
those  of  Fig.  5.5b.  Computations  of  p,  z  and  the  resulting  detection/false 
detection  probabilities  yielded  the  following  values: 

FIGURE  5.6  FIGURE  5.7 

(Aw  =  0.2  in  I  FOV  1 )  (Aw=0.2  in  IF0V2) 


Threshold  on  t 

0.4 

0.4 

new  M 

927 

871 

P  (D  ^  0) 

0.731 

0.857 

Z'NESR(D  i4  0) 

1 .81  x  lo"9 

2.33  x  10 

(D  =  2.66  x  1011; 

fSI 

II 

O 

NJ 

■tr 

for  10  second  dwell) 

p  (D  *  0) 

-0.288 

0.294 

z-NESR(O-O) 

1 .29  x 10'9 

1 .26  x 10 

P(TiT)  \  for 

~0.66 

'0.97 

P(T|0)  J  threshold 

0.075 

'0.05 

'0.07 

P(TlT) l  for 

'0.60 

'0.95 

P (T [ 0)  /  threshold 

0.08 

'0.035 

'0.055 
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WQVENUMBE*  ( CM-1  ) 

Figure  5-6  Same  as  Fig.  5-4a,  except  that  spectral  elements  corresponding 
to  t<0.4  have  been  assigned  a  contrast  value  equal  to  the 
mean  contrast  of  the  remaining  elements. 


*  ; 
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Figure  5.7  Same  as  Fig.  5.5a,  except  that  spectral  elements  corresponding 
to  t<0.4  have  been  assigned  a  contrast  value  equal  to  the 
mean  contrast  of  the  remaining  elements. 
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For  the  ideal  case  (no  interferences)  the  MDQ.  corresponded  to  Z  =  0.124 
and  R“0.075.  which  resulted  in  detection  and  false  detection  probabilities 
of  0.94  and  0.01,  respectively.  Thus,  it  is  interesting  that  only  slightly 
degraded  performance  is  obtained  in  the  non- idea  I  case  when  the  excess  H20 
and  the  emitting  target  S02  occupy  the  same  IFOV,  while  very  poor  perform¬ 
ance  is  obtained  when  they  occupy  different  IFOV's  (e.g.,  0.95/0.055  versus 
0.66/0.05).  It  is  interesting  because  it  is  really  the  former  situation 
that  led  us  to  include  spectral  masking  in  the  correlation  processing.  The 
explanation  of  these  results  is  that  weak  atmospheric  absorption  lines 
become  important  after  the  effects  of  the  strong  lines  have  been  eliminated. 
That  is,  weak  atmospheric  lines  within  the  S02  band  tend  to  add  to  the  $02 
lines  if  the  S02  is  observed  in  emission  and  the  excess  H20  is  in  IF0V2  (see 
Fig.  5- 3c),  while  they  tend  to  attenuate  the  S02  lines  if  the  excess  H20  is 
in  the  other  IFOV  (see  Fig,  5.3b).  In  fact,  the  opposite  of  this  "interesting" 
outcome  would  obtain  if  the  S02  band  is  observed  in  absorption.  In  any 
case,  it  is  evident  that  the  S02  amount  (MDQ)  must  be  increased  to  ensure 
that  adequate  performance  is  obtained  in  all  possible  situations. 

We  recalculated  the  detection  parameters  for  the  "Figure  5.6  case"  with 
D  increased  by  the  factor  1.5.  All  results  were  the  same,  except  p(D/0) 

—  --  ,  - Q 

which  increased  to  0.869,  and  z*NE$R(Dj<0)  which  increased  to  2.49x10  ; 

the  corresponding  P ( T ) T)  vs.  P (T | 0)  tradeoff  improved,  as  indicated  by  the 
following  set  of  values: 

R(threshold  on  [ p | )  P(T|T)  P (T 1 0) 


0.075 

'0.99 

-0.05 

0.08 

-0.985 

-0.035 

0.1 

-0.92 

0.005 
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These  values  suggest  that  the  MDQ  for  the  non- ideal  case  is  approximately 
1.5  times  that  for  the  ideal  case. 

We  also  recalculated  the  detection  parameters  for  the  "Figure  5.6  case" 
witn  the  additional  assumption  that  the  predicted  tgc  spectrum  is  based 
on  an  estimate  of  H^O  column  thickness  that  is  too  high  by  0.4  pr  cm.  Again, 

D  was  taken  as  1.5  times  the  original  MDQ.  The  resulting  P(T|T)  and  P(T|0) 
for  R* 0.075  are  '0.985  and  '0.055.  respectively.  Thus  we  can  say,  generally, 
that  real  world  effects  do  not  increase  the  MDQ  of  SOj  by  more  than  a  factor 
of  1.5.  if  spectral  masking  is  used  in  the  correlation  processing,  and  if 
pre-whitening  can  be  assumed  effective  in  essentially  eliminating  the  effects 
of  background  temperature/emi ss iv? ty  differences  in  the  two  IFOVs. 

5-6.2  The  HC1  2900~ 3000  cm  ^  Detection  Band 

The  normalized  contrast  spectral  radiance  ra+t/D  for  D  *  2.5  x  10^ 
(molec/cm  ) (W/cm  sr  cm  )  is  shown  in  Fig.  5.8a.  This  value  of  D  is  the 
MDQ  for  this  band  for  a  dwell  time  of  10  seconds.  Spectral  elements  were 
masked  out  if  their  predicted  transmittance  was  less  than  0.45,  which  reduced 
the  effective  number  of  elements  from  1000  to  617.  Again,  the  plotted  con¬ 
trast  of  the  masked  elements  is  equal  to  the  mean  contrast  of  the  remaining 
elements.  Of  course,  for  HC1,  which  has  a  very  sparse  emission  spectrum,  it 
would  have  been  better  to  use  the  HCl  spectral  absorption  coefficient  as  a 
criterion  for  masking,  but  it  was  easier  to  use  the  same  type  of  criterion 
that  was  applied  to  SOj.  Figure  5.8b  shows  the  normalized  reference  contrast 
spectrum  to  plotted  to  the  same  scale.  Figure  5.8  represents  the  case  of 
0.2  pr  cm  excess  HjO  in  IF0V1. 
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5.8  Part  (a):  Normalized  contrast  spectral  radiance  to  +  £/D  for 

.  1  \ 
the  HC1  2900-3000  cm  detection  band  for  0*2.9x10 

2  2  -1 

(molec/cm  ) (W/cm  sr  cm  ).  IF0V1  contains  0.2  pr  cm 
more  HjO  than  IF0V2. 

Part  (b) :  The  normalized  spectral  contrast  to  plotted  to  the 
same  scale  as  part  (a). 


Figure  5-9  shows  the  corresponding  normalized  contrast  spectra  when  the 
excess  HjO  is  in  IF0V2,  with  the  HC1,  rather  than  in  IF0V1.  The  number  of 
spectral  elements  that  survived  the  masking  criterion  is  606. 

The  detection  and  false  detection  probabilities  for  the  interference 

cases  represented  in  Figs.  5.8  and  5.9>  were  found  to  be  approximately  0.9*5/ 

0.005  and  0. 9*5/0. 007,  respectively,  for  a  threshold  value  R®0.1.  These 

values  may  be  compared  to  the  0.955/0.02  probabilities  obtained  in  the  ideal 

case  for  R*  0.075  with  no  spectral  masking.  This  does  not  imply  that  the 

interferences  Improve  performance,  but  rather  that  spectral  masking  would 

have  resulted  in  slightly  better  HC1  detection  performance  in  the  ideal  case 

2 

(as  a  result  of  increased  variance  Oq  for  the  reference  spectrum).  These 
results  show  that  interference  type  b  has  an  insignificant  effect  on  the 
MDQ  for  HC1. 

5-6.3  The  HBr  2450-2650  cm  *  Detection  Band 

Figure  5- 10a  gives  the  normalized  spectral  contrast  to+c/D  for 
HBr  when  the  quantity  present  is  the  MDQ  for  t^^lO  seconds,  and  the  excess 
H2O  is  in  IF0V1.  The  selected  transmittance  threshold  is  0.65,  which  reduces 
the  number  of  spectral  elements  from  2000  to  1 805 •  Figure  5.10b  gives  the 
reference  spectrum  tou 

For  the  ideal  case  with  no  spectral  masking  the  MDQ  corresponds  to 
P(T|T)  -0.9*2  and  P(T|0)-0.008  for  a  threshold  R-0.06.  The  detection  and 
false  detection  probabilities  corresponding  to  the  interference  case  repre¬ 
sented  in  Fig.  5.10  are  approximately  0.92  and  0.02.  Since  these  values  repre¬ 
sent  acceptable  performance,  we  can  say  that  interference  type  b  has  an 
insignificant  effect  on  the  MDQ  of  HBr  when  the  revised  correlation  proce¬ 
dure  is  used.  The  revised  method  is  effective  in  spite  of  the  fact  that  the 
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Figure  5.10  Part  (a):  Normalized  contrast  spectral  radiance  to  ♦  £/D  for 

the  HBr  2450*2650  cm  '  detection  band  for  0*2.44x10 
2  2  - 1 

(molec/cm  ) {W/cm  sr  cm  ).  IFOVI  contains  0.2  pr  cm 
more  1^0  than  IF0V2. 

Part  (b) :  The  normalized  spectral  contrast  to  plotted  to  the 
same  scale  as  part  (a). 
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transmittance  envelope  has  a  large  variation  over  the  blue  wing  of  the 
4.3  um  CO2  band  (2450-2530  cm  ')  (see  the  HBr  transmittance  plot  in  Vol.  II). 

It  should  be  noted  that  a  significant  part  of  the  interference  is  attri¬ 
butable  to  the  N20(2v^)  band  between  2525  and  2590  cm  This  band  appears 
in  the  radiance  contrast  spectrum  5  because  of  the  different  amounts  of  H20 
below  the  ^0  in  the  two  IFOV's.  However,  the  effect  of  the  N20/H20  is 
small  when  spectral  elements  are  excluded  on  the  basis  of  the  predicted 
transmittance.  For  both  the  HC1  and  HBr  detection  bands  a  modest  improvement 
in  performance  could  be  realized  by  excluding  from  the  correlation  nearly  all 
spectral  elements  not  containing  HC1  or  HBr  lines.  This  improvement  would  be 

realized  in  both  the  ideal  and  non-ideal  cases,  since  its  principal  effect  is 

2 

to  increase  the  apparent  variance  Og  of  the  reference  spectrum  ra.  The  same 
argument  applies  to  all  species  with  sparse  line  spectra. 

For  HBr  we  also  computed  the  effects  of  different  background  tempera¬ 
tures  in  the  two  IFOV's  (without  any  difference  in  the  H20  column  thicknesses) 
The  computation  is  intended  to  show  that  prewhitening  of  the  observed  contrast 
radiance  spectrum  is  an  important  part  of  the  spectral  correlation  processing. 
Figure  5.11  shows  the  normalized  spectral  contrast  to + £/D  when  the  back¬ 
ground  surface  in  IF0V1  is  I0°C  cooler  than  the  background  in  IF0V2  (290K  vs 
300K) .  The  same  spectral  elements  have  been  excluded  as  a  result  of  threshold 
ing  on  the  predicted  transmittance,  but  now  the  spectral  radiance  differences 
in  the  NjO  band,  in  the  H20  lines  above  2590  cm  ' ,  and  in  another  atmospheric 
band  below  2490  are  very  large  as  a  result  of  the  large  background  temperature 
difference.  There  is  also  a  large  variation  in  the  continuum  level  of  lj D 
that  results  from  both  the  variation  in  atmospheric  transmittance  envelope 
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Figure  5.11  Same  as  Fig.  5.10,  except  the  IFOV's  have  equal  column 

thicknesses  of  H.O  and  surface  temperatures  that  differ  by 
10  degrees  (290K  for  IFOVI  and  300K  for  IF0V2). 
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and  the  different  Planck  function  spectral  slopes  for  the  two  background 
temperatures.  Figure  5- 11b  shows  the  reference  rot  spectrum  plotted  to  the 
same  scale  as  Fig.  5.11a.  If  detection  of  HBr  were  based  on  the  correlation 
of  Figs.  5.11a  and  5.11b  (with  no  attempt  to  prewhiten  Fig.  5.11a)  the 
detection  and  false  detection  probabilities  would  be  approximately  0.66  and 
0.023,  respectively.  In  other  words,  the  occurrence  of  10-degree  surface 
temperature  variations  over  a  distance  of  10  meters  would  lead  to  unaccept¬ 
able  performance  unless  an  effective  prewhitening  scheme  is  implemented  prior 
to  spectral  correlation. 

The  general  features  of  Fig.  5* Ha  suggest  that  standard  methods  of 
prewhitening  would  not  be  effective  unless  they  are  preceded  by  somewhat 
elaborate  methods  of  spectral  rejection.  Effective  spectral  rejection  methods 
can  probably  be  devised  on  the  basis  of  a  priori  knowledge  (e.g.,  the  predic¬ 
tions  of  FASC0D1 ) .  It  is  also  probable  that  useful  spectral  rejection 
methods  must  be  devised  individually  for  each  of  the  detection  bands  of 
interest. 
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APPENDIX  A 


PROOF  OF  EQUATIONS  10  AND  U  IN  SECTION  2 

As  stated,  we  will  assume  the  atmospheric  spectral  transmittance  and 

gas  spectral  absorption  coefficient  can  be  predicted  exactly;  i.e.,  t'ol' 

v  gv 

can  be  replaced  by  T^a  in  Eq .  (9)  of  Section  2.  We  will  still  denote  the 

spectral  mean  of  t 'a'  =  t  ct  by  y',  and  the  mean  of  AN  by  y.  With  the 
v  gv  v  gv  v 

above  assumption,  y  =  Dy"  (see  Eq.  8  in  Section  2). 

Let  o'  denote  the  standard  deviation  of  the  reference  spectrum: 


o' 


[« 


t  a  -y 
v  gv  M 


Equation  (9)  of  Section  2  can  then  be  written 


(A-l) 


P 


<(ANv-y)(Tvagv-y')> 

r  9  1,/2 

0'[<(ANv-y)2>] 

<[D(Vgv~^)  +nv3(TvV^> 
a'{<[D(Tvagv-y')  +nv]2>} 


(A-2) 


where  the  second  form  results  from  substituting  Eq .  (8)  and  y  =  Dy'  into 
the  first  form. 

Expand  the  numerator  product  and  squared  denominator  term  in  (A-2), 
and  then  evaluate  the  spectra!  averages  term-by-term.  The  result  is  siiown 
in  (A-3)  where  we  have  used  (A-l)  and  written  the  spectral  averages  in 
terms  of  sums  over  M  resolved  spectral  elements  (the  particular  spectral 
elements  that  would  be  selected  in  the  band  of  a  given  species  to  provide 
optimal  detection). 


D<((t  a  -p')2)>+<\n  (t  a.  ~u')^> 
v  g\>  v  v  gy 

a'  D2<(t  a_  -y ')  2>+ 2D<n  (t  a  *u')^+<n  2>] 

L  N  v  gv  v  v  gv  '  '  v  'j 

a-[D2{a')2  +  20  J  nv(Tvagv-u')  +~  I  ny2] 

Equation  (A-3)  for  p  can  be  simplified  by  making  use  of  the  general 

rule  for  determining  the  statistics  of  a  weighted  sum  of  independent  normal 

variates:  If  v^  ,  v^,  vp  are  independent  normally  distributed  random 

variates  with  means  and  standard  deviations  (y,,a.),  (y_,a0),  ....  (y  ,o  ), 

it  L  L  n  n 

then  the  weighted  sum  v  =  ag  +  a  jV  j  +  ^-^2  +  •  •  •  +  anvn  '  s  norrnai  ly  d  istr  i  buted 
with  mean  and  standard  deviation 


m  *  a0  +  Vl  +  V2  +  +  Vn  * 

,  2  2  .  2  2  .  2  2  x1/2  ,,  n 

a  »  (3,  o,  +  a2  o,  +  ...  +  an  oR  )  .  (A-!*) 


Thus,  the  weighted  sum 

q  s  l  nv(Tvagv"»,') 


(A-5) 


appearing  in  the  numerator  and  denominator  of  (A~3)  is  a  normal  variate  with 
zero  mean  (since  the  n^  have  zero  means)  and  standard  deviation 


a 


[l  <\agV-v')2  an  J 


1/2 


0  H 

n 


1/2 


l  <Tvy^>2 


1/2 


rT  NESR  M,/2 


(A-6) 
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where  o  ■  fl  NESR  is  the  standard  deviation  of  the  noise  n  . 
n  v 

Now,  consider  the  scaled  sum  s  of  independent  unit  normal  variates 
x;  (which  have  zero  mean  and  unit  standard  deviation): 


_  M 

s  -  /2  NESR  o'  Ex,.  (A-7) 

i  -  1  1 


By  the  addition  rule  for  normal  variates,  s  is  normal  with  zero  mean  and 
standard  deviation 


o 


/2  NESR  o'  M1/2  ; 


(A-8) 


i.e. ,  the  random  variate  s  is  statistically  identical  to  q,  which  means 
that  the  sum  (A-7)  can  be  used  in  place  of  (A-5) .  We  note  also  that  any  of 
the  random  noise  samples  n^  is  Identical  to  NESR  x{ ,  where  is  unit 
normal;  thus,  the  last  term  in  the  denominator  radical  in  (A-3)  can  be 
written 


M 


l 

?  -  1 


x 


2 

i 


(A-9) 


Substituting  s  for  q  in  (A-3),  and  substituting  (A-9)  in  the  denominator, 
we  obtain 


P 


p(»-)2  +  £ 


M 


l 

-  1 


t>V)2  + 2D  *  ySR  <■'  I  K.  i  I 

M  i-1  '  M  i-l 


(A-10) 


133 


For  0  =  0  (no  target  gas  present) 


p(D-O) 


(A- 1! 


which  is  the  same  as  Eq.  (11)  of  Section  2. 

Note  that  the  bracketed  term  in  the  denominator  of  (A- 10)  is  the  sum 
of  a  perfect  square,  equal  to 


1 

M 


l  DZ(o')2  +  2  (Da')  /2  NESR  x.  +  2(NESR)2 
i  =  1  L 


- 

m 


i  M  r 

T. I, 

«  *  1  L 


Da'  +  Jl  NESR  x. 


2 

Using  this  fact  and  dividing  numerator  and  denominator  by  D(o')  ,  we  find 
that  for  0^0, 


(A- 13) 


where  Z  =  (a'D)/(/2"  NESR)  is  the  ratio  of  rms  signal  to  rms  noise  in  the 
contrast  radiance  spectrum  AN^.  This  result  is  the  same  as  Eq.  (10)  of 


13* 


Section  2. 


APPENDIX  B 


VALIDITY  OF  THE  APPROXIMATION  D*~  F# 


From  Eqs.  (22)  thru  (2k)  of  Section  2.5,  It  Is  seen  that  the  approxi¬ 
mation  requires  that 


(2F#)2  I  +  (2F^)2 


(B-l) 


for  some  choice  of  the  constant  C  (assuming  the  optics  and  scene  temperature 
are  equal).  If  we  let  t  *  0.261  and  e  ■  0.249,  and  choose  C  such  that  the 

u  h 

error  In  (B-l)  Is  zero  at  F  ■  l,  then  the  error  for  F  ' s  up  to  5.0  is  as 
shown  in  Table  B-l.  Note  that  the  maximum  error  is  five  percent. 


Table  B-l  Error  In  D*  Values  Scaled  Linearly  with  F* 


Percent  Error  In  D  ~  Fw 
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